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I. INTRODUCTION 



The idea that our 4-dimensional observed universe may be a hypersurface or "brane" in a higher dimensional 
spacetime is motivated by string- and M-theory |^. In particular, 5-dimensional braneworld scenarios, in which 



our universe represents the boundar" 



i§^a i,[i?,FTi^^,p,[i4iraF 

y, |4^, ^ 1^ ^, 



'EM 
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received considerable attention 




of a 5-dimensional s 



The case in which the bulk spacetime 

is Anti de Sitter space and orbifold compactification is realized with the brane as fixed point has been particularly 
studied. In this situation, it has been shown that 4-dimensional gravity is recovered on the brane pT| , pO| , |2^ , 23 on 
energy scales much lower than the brane tension and /or bulk curvature, and late time cosmology is not changed if 
the brane tension is sufficiently high |l^, 13, ^ |l^, |l^. In an attempt to solve the fine-tuning problem between 
the bulk cosmological constant and the brane tension, more complex models in which the bulk or the brane are filled 
with several species (such as scalar fields) have recently been proposed (see for example and references therein). 

In these models Z2 symmetry is often assumed, and this is particularly convenient when considering boundary 
conditions on the brane. If Z2 symmetry is dropped, brane motion in the bulk must be taken into account and 

mmetry 
there 



involved calculations are required in order to determine the boundary conditions on the brane. Whilst Z2 S" 
is motivated by M-theory and is required for a supersymmetric brane configuration, such as a BPS state 
exist situations in which Z2 symmetry is broken. This occurs, for example, when the brane is charged and c oup les to a 

^ , , . , , rim — I , .... rCm rm 1^ _j 1 — 1 1 1 rrrrim 1 — 1, 



4- form field in the bulk Cosmological asymmetric brane models have been studied in |2^, |2|, ^|3|, Hj. 

These developments have prompted us to derive gauge invariant perturbation theory for brane cosmology with one 
codimension. Our aim is develop a formalism which may then be applied to any situation of cosmological interest. 
Previously, perturbations in braneworld cosmology have been extensively studied in the literature mostly for the case 

H, |o[ mn, ^, n n, 0, 0, n, g, |oi ig. 



of Z2 symmetry [p5|, 
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Here we consider the most general 

situation in which the spatial background geometry on the brane has maximal symmetry and thus represents a space of 
constant curvature k. We do not assume Z2 symmetry, and the boundary conditions on the brane are discussed. Also, 
no particular gauge choice for the metric component g^^ is made. The perturbation equations in the bulk and on the 
brane are derived for general bulk and brane stress-energy tensors. This makes our formalism particularly convenient 
when analyzing situations in which different bulk components (such as several scalar fields) are also considered. The 
formalism can be used to study phenomena which have important observational consequences, the most important of 
them being the calculation of the anisotropics of the cosmic microwave background |l| . Since one must in general 
first determine the behaviour of perturbations in the bulk before being able to determine their behaviour on the 
brane |pj|, we pay particular attention to the relation between bulk and brane gauge invariant perturbation variables. 
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These become more subtle when the position of the brane is displaced. Indeed we define a set of gauge invariant 
variables in which the perturbation equations on the brane become similar to the usual 4-dimensional equations. We 
then study the new terms arising in braneworlds. 

Since we assume very general background spacetimes and no symmetry, some of our equations are extremely 
cumbersome. In order to guide the reader through the rest of the paper, we now give a general overview of the 
methods we use, the variables we introduce, and the equations we derive in this paper. 

The basic setup is one of a 3 + 1-dimensional brane where the 3-space of constant time is maximally symmetric 
(a space of constant curvature), embedded in a 4 + 1-dimensional bulk. As Zi symmetry is not assumed, the bulk 
spacetimes on each side of the brane will generally differ. Both the brane and the bulk may contain arbitrary matter. 
Our notation is as follows: 

• x", a = 0, 1, 2, 3, 4 : spacetime coordinates (Greek indices), with metric gaii and covariant derivative Dq., 

• x% I = 1, 2, 3 : coordinates on the maximally symmetric 3-space (second part of Latin alphabet) with metric 
7ij and covariant derivative V;, 

• ct", a = 0, 1, 2, 3 : brane-worldsheet coordinates (first part of Latin alphabet), 

• X"((T°) : brane position in target-space. 

• A Roman subscript b indicates "brane" whilst B denotes "bulk" . 

• Certain variables such as the brane matter content (/^(cr"), p((t''), etc) are only defined on the brane. Other 
variables such as the normal vector to the brane _L" or the extrinsic curvature /Ca/3 are also defined at the brane 
position, but since they describe the embedding of the brane in the bulk, they may take different values on 
either side of the brane. All these brane-related variables are underlined. 

The action for the system is 

S = Seh + S'^ + + Sqh 

d'x^l (^^i? + C"^^ + J d'a^lC^ + Sgh. (1.1) 

Here Squ is the Gibbons-Hawking boundary term required to consistently derive the Israel junction conditions [ p4| , 
and K5 is the fundamental 5-dimensional Newton constant (related to the 5-dimensional Planck mass M5 by K5 — 
Gtt^MI). Furthermore, R is the bulk scalar curvature, gap and Xab are the bulk metric and the induced metric on 
the brane respectively, and and £™ are respectively the Lagrangians for arbitrary matter in the bulk and matter 
confined on the brane. They may also contain a cosmological constant or brane tension. The induced metric on the 
brane is (see for example ^T\) 

Xafc(a) = <?^.(X)^^, (1.2) 

and the Einstein equations resulting from action ( |l.l| ) are 

Gaji = K5 (T'q/s -I- DT^^) , (1.3) 



where D is a covariant Dirac (5- function specifying the position of the brane (see Section [II B ) , and 

~d^~da^ \Sxab 



^^(S,], (1.4) 
r^^XsSrf^). (1.5) 



As noted above, we underline T_"^ and Xab t o emphasize that they are only defined on the brane (see Section HI). 

We consider these Einstein equations ( |l.3| ) for a homogeneous and isotropic brane and bulk background with first 
order perturbations. As is summarised schematically in the left hand panels of Fig. |l|, these equations contain three 
parts: one is continuous; the second is discontinuous across the brane; and the third part is singular at the brane 
position (proportional to D). The coefficients of each of the individual parts must be equated, leading to a number of 
different equations. The continuous part gives the Einstein equations in the bulk and, via the Gauss-Codacci equation, 
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they also determine the 4-dimensional Einstein tensor on the brane (see Section LV). The discontinuous (but non 
singular) part is only non-trivial when Z2 symmetry is not assumed. It then gives equations for the continuous part of 
the extrinsic curvature, and it describes the energy and momentum exchange between the brane and bulk, leading to 
the equation of motion for the brane — the so-called "sail equation" |T[ |3^ . Finally, the singular part represents 
the second junction condition which relates the bulk geometries on each side of the brane through the brane geometry 
and matter content. 

We also discuss the so-called conservation equations for the stress-energy tensors given in Eqns (1.4,L5|). Again, 
these contain a discontinuous, continuous and singular part. As is summarized schematically in the right hand panels 
of Fig. 1^, the continuous part gives the bulk energy momentum conservation, the discontinuity simply describes the 
conservation of the jump of the bulk stress-energy on the brane, and the singular part leads to energy momentum 
conservation of the brane with a possible contribution from the bulk, and to the sail equation. 

When discussing the perturbations of these equations, we will make use of the maximal symmetry of the 3- 
dimensional subspaces parallel to the brane. Our geometrical quantities will be decomposed into scalar, vector and 
tensor degrees of freedom (with respect to these 3-spaces). This decomposition is not identical to the more physical 
one containing density modes, vortic ity m odes, and 5-dimensional gravitational waves. The relationship between these 
two approaches is given in Section VA. Finally, in order to set up a consistent gauge-invariant formalism for the 
evolution of these perturbations, we will see that it is crucial to take fully into account the perturbed brane motion 
(which can be written in a gauge invariant manner). This degree of freedom will be central to our analysis. 

The outline of the paper is the following. In the next section (Section ||) we discuss the unperturbed (or background) 
5-dimensional bulk: we allow a foliation (with two codimensions) into maximally symmetric 3-spaces, and do not 
specify the presence of the brane. The Einstein and conservation equations for the bulk background are derived. In 
Section [II we introduce the brane and we discuss the boundary conditions at the brane position for the unperturbed 
spacetime without imposing Zi symmetry. In Section IV we derive the background equations for an observer on the 
brane. In Section we perturb the background. We introduce gauge invariant variables and derive the perturbed 
Einstein and conservation equations in terms of these va riab les. The perturbed brane including the perturbation of 
the brane position is discussed in Section |v^ . In Section VII we reformulate the perturbation theory from the point 
of view of an observer confined to the brane, and in the last section we draw some conclusions. 

Finally, we also provide an extensive and highly technical Appendix where we present all the relevant intermediate 
steps required to obtain the results presented in the text. (Examples are, for instance, the perturbed Christoffel 
symbols and the components of the perturbed Riemann and Weyl tensors.) The Appendix is, in fact, more general 
than the main text since there we consider an + 1-dimensional brane (with an A^-dimensional maximally symmetric 
subspace) embedded in a A^ -|- 2-dimensional bulk: in the text we have set A^ = 3. Furthermore, whilst the text 
presents the perturbation equations in full generality, some specific examples such as a bulk scalar field are discussed 
briefly in the Appendix. 



II. BULK BACKGROUND 



In this section we describe the bulk background geometry and energy content without introducing a brane. We 
assume that the space orthogonal to the fifth dimension is maximally symmetric so that a homogeneous and isotropic 
brane can be accommodated, as discussed in the next section. We consider the most general stress-energy tensor 
which satisfies these symmetry conditions, and then derive the Einstein equations and the conservation equations. 



A. Metric and notation 



We consider a 5-dimensional spacetime with one timelike coordinate = ry, and four spacelike coordinates 
{x^, x^, a;'^, x**}, where = y. We assume that the constant time hypersurfaces are locally of the form x R, 
where is a 3-dimensional maximally symmetric space, i.e., a 3-space of constant curvature, parameterized by 
the coordinates {x^,x'^,x^}, with spatial metric a?{ri,y)^ij. The curvature of this space will be denoted by k. For 
example, we may choose the coordinates {x^, x^, a:'^} such that 

kx^X"^ 

where 5ij is the Kronecker symbol. The last spacelike coordinate y (the "extra dimension") is orthogonal to the 
maximally symmetric space. The metric has the signature H . The line element of the metric can therefore 
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FIG. 1: Structure of the Einstein equations and of the energy momentum conservation equations in the coordinate system (2.2), 
where coordinates 0, i are also brane coordinates and 4 represents the direction orthogonal to the brane. In components, the 
Einstein equations can be split into three parts: {fJ^v}, {/i4}, and {44}, where /x, v run on indices 0, i. These three parts possess 
a continuous part (defined everywhere in the bulk) and a jump at the brane position. Part {fiv} also exhibits a singular term 
at the brane position. The role played by all these terms is shown in the above diagrams. 
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be written as 

ds^ = n^dr/^ - a\jdx'dx^ - ft^dy^ (2.2) 

An overdot will denote derivation with respect to t], and a prime derivation with respect to y. In addition, we shall 
define 

du = -dr,, (2.3) 
n 

dn = \dy. (2.4) 

Covariant derivatives with respect to the full metric will be denoted by Da., and those with respect to 7jj by Vj. For 
convenience we also define 

n^^-'- , X. 1^ , (2.5) 

na nn no 

la' In' lb' , , 

H=-- , I = , U=--, 2.6 

a on 



as well as 

V' = ViV' , , Vi,=ViV,- , K=\. (2.7) 

Notice that du and dn do not commute since dudn — dndu = Idu — Udn- 

B. Stress-energy tensor 

We now consider the bulk stress-energy tensor T^p whose energy flux need not be at rest with respect to our (ry, y) 
coordinates. Let C/" be the normalized timelike eigenvector of with eigenvalue po- Correspondingly let N°' be the 
normalized eigenvector orthogonal to both and to the maximally symmetric 3-spaces, with eigenvalue Yq. Finally, 
let Pq be the eigenvalue of the three eigenvectors parallel to the maximally symmetric 3-dimensional slices. Note that 
any symmetric tensor can be decomposed in this way, and that the symmetry requires that the eigenvectors parallel 
to the symmetric 3-spaces are degenerate. In terms of these variables, the bulk stress-energy tensor can be written as 

T„/3 = {Po+Po)UJJi3-{Po-Yo)N^Np-Pogai3, (2.8) 
where the vectors C/" and N"^ are given by 

= (^'^'0'^^^) ' U^U>' = 1, (2.9) 
iV" = (-^/^^'O'-^^) ' N^N^' = -l. (2.10) 

Here /3 represents the Lorentz boost which must be performed along the y axis in order to be in the fluid's rest frame. 
When one is not in the rest frame of the fluid, both its energy density and its pressure along the extra dimension are 
modified, and the fluid exhibits a flux through an ?/ = constant hypersurface. As usual, 7 = l/-\/l — 0"^. 

Below it will be more convenient to use a different definition for the stress-energy tensor componcnits a definition 
which is less adapted to the fluid, but better adapted to our coordinates. To derive it, let us denote by the 
5-velocity of a bulk observer who is at rest with respect to our coordinate system, 

««=(^io,0^ , uX = l' (2-11) 

and by the spacelike unit vector orthogonal to both u" and M., 

0,0, , n^n'' = -l , nX = 0. (2.12) 
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Note that neither nor n" are geodesic vector fields, but their orthogonality and normalization is conserved. The 
most general form of the bulk stress-energy tensor satisfying the required symmetry with respect to translations and 
rotations in can be written as 

Tai3 ^ {P + p)uaU(} - {P - Y)nanp - Pgap ~ 2Fu(^anp), (2.13) 

where f(a9i3) = \{fagp + gafp) denotes symmetrization. In components this gives 

Too = n'p, (2.14) 

T^j = a^Pj^j, (2.15) 

ro4 = -nbF, (2.16) 

T44 = b^Y. (2.17) 

Thus p — T^^u^u^ is the bulk energy density as measured by an observer with 5-velocity , F = T^^u^'ny is the 
energy flux transverse to M.^ and P, Y are the pressure along the directions a;*, y, respectively. The new variables p, 
Y , P and F are related to the old ones by 

P = 7'(Po+/3'Fo), (2.18) 

Y = 7'(>o+/3'po), (2.19) 

P = Po, (2.20) 

F = /37'(po + lo). (2.21) 

The last relation again shows that F represents the energy flux in y direction. This flux, as well as the energy density 
p and pressure Y along the y direction measured by an observer at rest with respect to the coordinate system, are 
obtained from the components of the stress-energy tensor in a frame at rest with respect to the fluid simply by a 
Lorentz transformation. Somewhat more complicated but equ ally straightforward expressions express the old variables 



in terms of the new ones (see Appendix C2). Note that in (2.9) the bulk velocity of the fluid is /? = n^JJ^ /u^U" 
When /3 = 0, 7 = 1, we recover the case in which [/" = and A^" — n", so that the rest frame of the bulk matter 
and the coordinate system coincide. 



C. Einstein equations 



Th e Christoffe l sym bol s, th e Riemann, Ricci, Einstein and Weyl tensors for the metric (2.2) are given in Appen- 
dices B3, B4, B6, 87 and B8 respectively, and the background bulk Einstein equations are 



With the stress-energy tensor ( 2.13 ) and the Einstein tensor from Appendix B6, Eq. ( p.22| ) becomes 

iK + m{'H+U) -i {dn + 2H)H = H5P {00} , 

-K-3{H^-H^) - [du +U){U + 2n) + {dn + /)(/ + 2H) = n^P {i]} , 

3 [duH + nH- HI) = K5F {04} , 

- 3A' -3{du + 2n)n + 3iH + I)H = n^Y {44} , 



(2.22) 



(2.23) 
(2.24) 
(2.25) 
(2.26) 



where we have indicated in braces on the right hand side from which component of the Einstein tensor these bulk 
Einstein equations are derived. Equations ( 2.23 , 2.26 ) were first discussed in [Q, and integrated with respect to the 
fifth dimension in [|l^, for the case of a negati ve bulk c osmological constant, P — Y ^ —p = A. 

The first and the third of these equations ( 2.23 2.25 ) are constraints (i.e., they do not involve second derivatives 
with respect to time). The other two are dynamical equations. In fact, there are only two independent dynamical 
variables which can be written as a combination of the scale factors n, a, and h. One can choose coordinates to remove 
this ambiguity: for example, in Gaussian coordinates & = 1 as in p^ , and in conformal coordinates b = n p9| , 
Of course other choices of coordinates are also allowed. We shall, however, keep b undetermined, so that any useful 
choice for b can be made at the end. 
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D. Conservation equations 

The Bianchi identities lead to the so-called conservation equations for the stress-energy tensor, 

JJ^J.^,a ^ Q (-2 27) 

Only for a = and a = 4 are there non-trivial relations, 

duP + m{P + p)+U{Y + p) + (dn + 3H + 2I)F = {0}, (2.28) 

{du + 3H + 2U)F + dnY + 3H{Y - P) + I{Y + p) = {4} . (2.29) 

These are the conservation equations for the energy density and the energy flux of the bulk components, respectively. 



The generalisation to several components is straightforward (see Appendix C4 ) . Written in term of the intrinsic fluid 



quantities, they give an equation of evolution for the energy density po and for fluid bulk velocity p. 

III. BULK BACKGROUND WITH A BRANE 

We now consider a homogeneous and isotropic 3-brane orthogonal to y (lying in the space of maximal symmetry) 
as a singular source, with intrinsic stress-energy tensor T^,^. 

A. Brane position, induced metric and first fundamental form 

Let us choose the intrinsic brane coordinates (cr°, ct') = (ry, x*), and embed the brane according to 

X" = r,, (3.1) 

X' ^ x\ (3.2) 

X'^ = yb = constant. (3-3) 

Note that it is always possible to choose the background coordinate y such that the unperturbed brane is at rest: this 
is the only coordinate choice made in this paper. 

As we shall see, the presence of the brane will introduce discontinuities at y = yb in several variables. For that 
reason, it is useful to decompose a given function / as 

/ = [/] {0{y - yb) - i) + (/) {y), (3.4) 

where 9 is the Heaviside function. This equation defines the continuous function (/) (y), whilst the discontinuity or 
jump of / when going from one side to the other side of the brane is given by 

[/] = lim (/(yb + £) - /(yb - e)) = /+-/". (3.5) 
Notice that we have the two product relations 

= {f){9) + \[f][9], (3.6) 

[fg] - (/>[.9] + [/](.9)- (3.7) 

For later convenience, and when considering a continuous function (/), we will also define the continuous part and 
the jump of its derivative by 

(9„)(/) ^ (a„(/)), (3.8) 

[9„](/) ^ [9„(/)]. (3.9) 

Sometimes we shall also need [/] for variables / describing the embedding of the brane, and thus which may take 
different values, /^, /^, on either side of the brane. The quantities [/] and (/) are defined by 

[/] = /+-/-, (3.10) 

(/) ^ k{t + r)- (3-11) 
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The normal unit vector to the brane, _Lq, is given by 



, = (3.12) 



One obtains (up to an overaU sign) 

(3.13) 

As we shall see, b can be discontinuous on the brane and _L" can have different values on either side of the brane. 

From the induced metric one can define the first fundamental form [ ]69| qa/s = gap + J^qJ^/j, where gap is again 
evaluated on (either side of) the brane, and we have QafiAJ^ — 0. On the brane, Qap is related to Xabi^') by 

We can decompose the stress-energy tensor on the brane, T^^(X), as 

Tap = (£ + P)llaUp - Pgap, (3.15) 
where is the 4- vector of the energy flux on the brane matter, 

u"=Q,O,0^ , u^u^ = 1. (3.16) 

Note that T_^^1J^ = Ik^X.^ = 0. This is the most generic stress-energy tensor compatible with a homogeneous and 
isotropic brane. 



B. Einstein equation 

In the presence of the brane, the 5-dimensional Einstein equations become 

Gap = K5 {Tap + ^TLap) > 



where, from Eqn (1.1), the "covariant Dirac function" D is 



D = 



5{y-yh)- 



(3.17) 



(3.18) 



Here g and q are the determinants of the metric gap and first fundamental form q ap re spectively, evaluated at the 
brane position. Written in components, the Einstein equations with the brane, Eq. ( |3.17D , become 



?,K + 3n{n+U) -3 {dn + 2H)H = 
K-3{n^- - {du +U){U + 2H) + {dn + /)(/ + 2H) = 

3{duH + HH - HI) = 
- 3K - 3 {du + 2n)H + 3 {H + I) H = 



K5{p + r>p) {00}, (3.19) 

^5 {P + DP) {ij} , (3.20) 

{04} , (3.21) 

K5Y {44} . (3.22) 



A global solution to these equations has been derived in |Tj| with the assumption of a pure ne gative cosmological 
constant in the bulk, and using Gaussian coordinates. The right hand sides of Eqns (3.19, ^^ ) contain a singular 
term proportional to D due to the presence of the brane. As we will see below, although the first fundamental form is 
continuous on the brane, its first derivative with respect to the fifth dimension y (i.e., the terms H and /) may jump 
and its second derivative {dnH and 9„J) can be singular. Thus the Einstein tensor contains a singular part which 
must be matched with the singular part of the stress-energy tensor. We now turn to the problem of relating these 
terms to the brane matter content. 



Forbidden 



Allowed 
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FIG. 2; Schematic illustration of the first Israel condition. We have embedded in a Minkowskian space a 2-dimensional spacelike 
bulk of metric ds^ — a^dx^ + b^dy^. The brane is the thick horizontal line in the middle of both panels, and the grids represent 
the metric coefficients so that the grid spacing is proportional to b and a along the vertical and the horizontal directions 
respectively. In the left panel, [a] 7^ 0, [b] — 0, in the right one, [6] 7^ 0, [a] — 0. The first Israel condition states that when 
considering a line of constant x, there must not be any discontinuity when crossing the brane: this is obviously not the case in 
the left panel. On the contrary, nothing is said about how the spacing of the horizontal lines evolves across the brane. This 
translates into the fact that 6 is allowed to be discontinuous (right panel). 



C. Israel junction conditions 

The extrinsic curvature formalism is a useful tool in the analysis of junction conditions on a singular surface [[70| . 
The first Israel condition [ ]65[ imposes the continuity of the first fundamental form, 

= 0. (3.23) 

Hence qa/3 is well-defined on the brane. Since we have goo = "^^(i'^) ^^'^ Iv ~ ~'3;^(X)7y, this condition implies the 
continuity of the scale factors n and a: [n] = [a] = (see Figure B). Note that the continuity of the metric function b 
is not required by the junction conditions and will not be assumed in what follows^ (see also Appendix^. 

Nevertheless, the first derivative with respect to y of a and n (which are proportional to / and H), are allowed to 
jump. In order to study the behaviour of these quantities on the brane we consider the extrinsic curvature tensor (or 
second fundamental form) with respect to the brane, namely 

!Cap=qlD^±^). (3.24) 

For the background metric, the components of the extrinsic curvature are 

£00 = -n^I, (3.25) 
ICij = a^Hjij. (3.26) 

Let us define the surface "stress tensor" S_af3 on the brane by 

Sa0 = T_al3 ~ -^TgafS- (3.27) 

Then the second Israel condition [6^ relates the jump in the extrinsic curvature with the energy content on the brane 
and requires that 

[Kap] = —i^bSxip- (3.28) 



^ Allowing b to be discontinuous makes the covariant Dirac function D ill-defined. This is not a serious problem, as all the terms involving 
this function can be grouped together to give the second Israel condition. Therefore, we shall continue to use the notation D and suppose 
that when b is not continuous, it corresponds to a regularized and mathematically consistent expression. 
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Allowed 




FIG. 3: Illustration of the second Israel condition. With the same conventions as in fig. ^ we show an example where [H] 7^ 0. 
This is possible if the energy density on the brane is non zero (it is positive in this illustration). 



(Note that the choice of the sign here is consiste nt w ith our choice for the sign oi in Eq. (3.13).) For our background 
this condition can be written as (see Appendix D7) 



[H] ^ -K5-p. 



p 



(3.29) 
(3.30) 



(See also Figure Bl) Alternatively, Eqns (3.29,3.30) can be obtained directly from the singular part of the Einstein 
equations (Is^IEM). 



D. Boundary conditions in the bulk 



We now comment briefly on the question of boundary conditions at the brane. Consider first the bulk Einstein 
). They form a system of second order partial differential equations in rj and y, and in order to 
solve them we must specify initial conditions on a spacelike Cauchy hypersurface, boundary conditions far from our 
braneworld (at infinity in a one brane scenario, or on another brane), and boundary conditions at our brane. The 
Israel junction conditions impose the continuity of a and n, and fix the jump in their normal derivatives at the brane. 
Since the Einstein equations represent a set of second order partial differential equations, these junction conditions 
are sufficient to allow us to solve the Einstein equations everywhere in the bulk. 
We now turn to the Einstein equations on the brane. 



equations (ET9-|!22 



IV. THE BRANE POINT OF VIEW 



An observer on the brane will not see 4-dimensional Einstein gravity. This may, however, be recovered in particular 
situations at low energy. The 4-dimensional Einstein tensor in general depends on bulk quantities and is quadratic in 
the brane stress-energy tensor. 

Here we discuss the 4-dimensional "Einstein equations" which lead to the modified Friedmann equations, and also 
the conservation equation on the brane. Finally we interpret the deviation from the 4-dimensional theory in terms of 
the 5-dimensional one. 



A. Einstein gravity on the brane 

Through the Gauss-Codacci equations, we can write the 4-dimensional Einstein tensor '■''•'G^^ in terms of the bulk 
stress-energy tensor T^p, the extrinsic curvature ICap and the projected Weyl tensor £_ap- The details of the calculation 
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can be found in ||3J], and the result is 

The projected Weyl tensor f ^.^ on the brane is obtained from the bulk Weyl tensor Cap^s as follows 



2 1 



(4.1) 



(4.2) 
(4.3) 



Here f[agp] = ^(/a5/3 — 9a fp) denotes antisymmetrization. This projection represents the contribution of the free 
gravity in the bulk to the gravity on the brane. In components, we have 



F 



—n Z, 
2 

1 2 



Z = K+{du+U){U-H)-{d„+I){I-H). 



(4.4) 

(4.5) 
(4.6) 



There is only one independent component in the Weyl tensor (as well as in its projection on the brane). This is 
related to the fact that this spacetime is the 5-dimensional analog of a 4-dimensional type- 1? spacetime in Petrov's 
classification IpSl. 



B. Friedmann equations on the brane 



Since the tensor ^'^■'G^^ contains only derivatives of the continuous first fundament al fo rm with respect to 77 and x% 
it is cont inuo us. Hence, on taking the continuous part of the right hand side of Eq. (4.1) (and applying the product 
relation (|3.6|)), we find the projected 4-dimensional Einstein equation on the brane. 



(4)^ 



.Q/3 



«5 (T^.)9^i£^ 



{T,.±^±n + 4 ) gap 



HQ {iCap) + im {K^p) + \qapmf - {Kn {k,.)) 



{£ap)^ 



(4.7) 



The right hand side of equation (4.7) can be split into four parts. The first depends on the average of the bulk stress- 
energy tensor, (Tap). The second is given by four terms quadratic in the average of the extrinsic curvature {JCap). 
These terms are known once the bulk Einstein equations have been solved. They vanish when Z2 symmetry is assumed 
— we return to this point below. Then there is a third part which contains four terms quadratic in t he jump of the 
extrinsic curvature Kq/j]. We have already determined these through the second junction condition (3.29 3.3C), and 
they are related to the brane stress-energy tensor: these terms will be responsible for the non-standard p"^ contribution 
in the brane Friedmann equations. Finally there is a fourth part, the average of the projected bulk Weyl tensor on 
the brane, describing the effect from the free gravity in the bulk. 

In components, we obtain the modified Friedmann equations which contain a dynamical equation and a constraint: 



-3{Hy 



(Z) 



- 2dun -m'-K = (p + p + 3y) + -^{p + 2P)p - {h) (2/ + h) + - (z) 

6 12 - — 6 

where we have isolated the continuous part of the projection of the bulk Weyl tensor, 

{z)^K-duH + {idu + u-n)u)-{dn + i) {i - h) - A (p + p). 



(4.8) 
(4.9) 

(4.10) 
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(Thes e equation s could alternatively have been obtained from the continuous part of the Einstein equations, 
Eqns (|]l9|-|^).) 



The cosniological consequences of these equations have been studied in [y_3|, |15|, 16, |17|, |1^ with assumption of Z2 
symmetry, in which case {JCap) — {H) — {I) = 0. These authors considered a negative cosmological constant in the 
bulk and assumed that the brane stress-energy tensor consists of a rigid part — the brane tension — and a fluid, 



T 



Ma0 + T 



f 



(4.11) 



In this case the bulk stress-energy tensor can be tuned to the brane tension in such a way that deviations from 
standard Friedmann equations are effective only at energies of or der A and h igher. 

If we assume F = (H) =Y — 0, integrate the sum of Eq. ( 3.21 ) and ( 3.22 ) once with respect to time, and compare 
the result with ( |4.8| ), one obtains an expression for the sum of the continuous part of the projected bulk Weyl tensor 
and the continuous part of the bulk energy density and pressure: this behaves as a radiation term. 



(Z) 



P)=C/c 



(4.12) 



where C is an integration constant (see Appendix E2) 
so that (Z) oc a^'* 



In P, |§ this is discussed for the case p + P^O 



Furthermore, notice that Eq. (4.11) is not a necessary requirement in order for the 5-dimensional Friedmann 
equations on the brane to reduce to the standard 4-dimensional Friedmann equations at low energy. It sufRces that 
the brane stress-energy tensor is dominated by a term which is almost a cosmological constant today (i.e., it can be 
a slow-rolling scalar field, see for example 123]). In this case, one has 



P = P4 

-Prk — A, 



(4.13) 

(4.14) 
(4.15) 



where pf , P/ are the energy density and pressure of the ordinary matter content on the brane. Now, if in addition 



\P + P\ « \Yl 
Y ^ A, 

A ~ iK5A^ 
6 



(4.16) 
(4.17) 

(4.18) 



the first two terms of the right hand side of the above Friedmann equations on the brane are proportional to K4p-^ , 
KiP_^ , with 



K4 = T-^S-^ ^ST' 



A 

A^ 



(4.19) 



and one re-obtains a te rm li near in the brane matter content. Note also that there are no extra A^ appearing in the 
projected Weyl tensor (4.10) since with these conditions P_ + pis A- independent. This was also noted in Ref. |]5^. 



C. Closing the system when Z2 symmetry is broken 



When solving the Einstein equations on the brane (4.7), we need the continuous part of the extrinsic curvature 



ilCaf})- If Z2 symmetry is assumed — as motivated by M-theory — the evolution of the bulk is the same on both 
the sides of the brane. In this case the Israel conditions determine £a/3 entirely: it is always possible to choose a 
coordinate system in which ri(yb + y) = '^^(yb — y) and a(yb + y) = a(yb — y), i-e., where n and a are even. Thus / 
and H are odd and the continuous value of lCai3 across the brane vanishes, 

{Kcp) - (Z2 symmetry). (4.20) 

This implies 

(/) = (F) = 0, {Z2 symmetry), 

= -I- = \ [I] > 

H+ = -H- = ^[H]. (4.21) 
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If Z2 symmetry is not assumed, as in tliis paper, the evolution on either side will in general be different and {ICafj) 
no longer vanishes. 



{Kcfi) ^ 



(Z2 symmetry broken). 



(4.22) 



One can, however, obtain a condition for the continuous part of the extrinsic curvature by considering the jump of 
Eq. (O). We obtain 



= 



(Notice that [^'^^G^p] 



This becomes, in components. 



[Sc. 



m{H)-p{i) = \ 



[p- 
[p- 



/3J • 

p + ZY] 



(4.23) 

(4.24) 
(4.25) 



where, using Eq. (4.6) and the junction conditions ( |3.29D and ( |3.3Cl| ), the jump of Z on the brane can be expressed as 



[Z] = {du + 2 (U) - H) [U] - [dn] {I - H) 



1^5 {I){2P+-p 



^5 {H)[P + -p 



(4.26) 



(Note that Eqns ( 4.24,4.25) could alternatively have been obtained from the discontinuous part of the Einstein 
equations ( |3.19|-|3.22|).) 

Equations ( |4.24 . 4.25D allow one to fix the unknown quantities (i/), (/), provided the jumps of both the bulk matter 
content and the Weyl tensor are known. Thus the continuous part of the extrinsic curvature depends not only on the 
brane matter content b ut also on the discontinuity of the bulk stress-energy and the projected Weyl tensors. If both 
vanish, Eqns (4.24 4.25) allow in p articular the trivial solution {H) = {I) — 0, which holds with Z2 symmetry. The 
jump of the Weyl tensor, Eq. ( 4.26 ), contains first derivatives of the extrinsic curvature with respect to y, and so it is 
not possible in general to determine (/) and {H) without first solving the Einstein equations in the bulk. Nonetheless, 
if the projection (4.26) of the 5-dimensional bulk Weyl tensor on the brane is known a priori (as in the case of a 
Schwarzschild- Anti dc Sitter bu lk with a known black hole mass on both side of the brane) , then (/) and {H) can be 
determined directly from Eqns ( 4.24 4.25| ) (see |^, ^ for a more detailed discussion). 



D. Brane motion 



On contracting Eq. ( [4.23| ) with the first fundamental form g"^ and using the second junction condition ( 3.2^ ), one 
obtains 

Z"" (iC„.) - U^^.T'^n . (4.27) 

This equation is known as the "sail equation" ^ The right hand side is an external force density on the 
brane due to the asymmetry of the bulk stress-energy tensor on the two sides. In analogy with Newton's second law 
(here the force is due to a pressure difference between the two sides of th e brane), T"'^, (JCafs), and [J^^J^i/T^''] play 
the role of mass, acceleration, and force, respectively. Notice from (4.20) that when Z2 symmetry is assumed, this 
equation vanishes identically. When Z2 symmetry is broken, the "acceleration" {JC^S) is non-zero. In this paper we do 
not assume Z2 symmetry, but recall that we have chosen a coordinate system in which the brane is at rest: Eq. (4.2J) 



must therefore be understood as dictating the condition that must be satisfied by {H) and (/) (and therefore b y the 
coordinate system itself) for the brane to remain at a fixed position yb. Later, however, we will see that Eq. (4.27) 
does indeed give a more intuitive equation of motion for the perturbed brane position or brane displacement (see 
Section ^rg). 



In components the sail equation leads to 



{I)p + -i{H)P=[Y] 



(4.28) 



which can also be obtained by taking the disco ntinuous p art of the {44} component of the Einstein equ ation. Eg . ( 3.22| ) 



or, of course, by a linear combination of Eqns (4.24 
does not involve the Weyl tensor 



4.25| ). This is the only combination of Equations (4.24,4.25) which 
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E. Conservation equations 



The singular part of the 5-dimensional energy conservation equation (2.28) yields the stress-energy conservation 
equation on the brane: we find 

duP + m{P + p) = -[F]. (4.29) 

(Again the generalisation to several interacting components may be found in Appendix D5.) Notice that the jump 
in the bulk energy flux transverse to the brane enters in the conservation equation, meaning that the brane matter 
content can act as a source or a sink to the energy flux along the flfth dimension. When this energy flux is continuous, 
the conservation equation on the brane reduces to the usual one, as discussed in [p4l . Another conservation equation 



appears in brane cosmology: by considering the singular part of Eq. (2.29) we obtain again the sail equation (4.28) 



Both equations were first found in |13] for the case a bulk cosmological constant. 



V. BULK PERTURBATIONS 



We now turn to perturbed quantities, and begin in this section by analysing the properties of the perturbed bulk: 
the perturbed brane itself will be introduced in Section We work throughout with gauge independent perturbation 
variables, which are inspired from a generalisation of the Newtonian (or longitudinal) gauge to the 5-dimensional case. 
First we introduce the bulk metric perturbation variables using the standard scalar, vector, tensor decomposition. 
We study their gauge transformation properties and define gauge invariant cornbinations. Then, in Section |VD , the 
perturbations of the bulk stress-energy tensor are considered, leading, in Section VE , to the gauge invariant perturbed 
bulk Einstein equations. Finally we write down the perturbed conservation equations (Bianchi identities). 



A. Classification of the perturbations 



Let us consider the perturbations of a spacetime with one timelike and n spacelike coordinates. The perturbed 
metric of this spacetime possesses ^(n -I- l)(n -I- 2) different components. Amongst these, a coordinate transformation 
allows n -|- 1 of them to be fixed, so that there are ^n{n + 1) independent metric coefficients. For example, in 
synchronous gauge, the Sgoa are set to zero. 

When solving perturbation equations about a given spacetime, one is naturally led to classify perturbations. Two 
classifications are of particular relevance. Firstly, the perturbations may be classified according to their physical 
meaning, and this is done by looking at the spin of the perturbation in a locally Minkowskian frame. The different 
perturbations are density (spin 0) modes, vorticity (spin 1) modes, and gravitational (spin 2) waves. Secondly, the 
perturbations may be classified more geometrically in terms of irreducible components under the group of isometrics 
of the unperturbed spacetime. This leads to scalar, vector and tensor perturbations. Under some circumstances, 
these two classifications are identical. In particular, this is true for a Friedmann-Lemaitrc-Robertson- Walker (FLRW) 
spacetime, which can be foliated by a set of maximally symmetric spacelike hypersurfaces. In brane cosmology, 
however, the bulk is not as symmetric as in the FLRW case, and the two classifications are different. 

Components which transform irreducibly under symmetries of the background spacetime evolve independently (to 
linear order) while the physical spin components mix. 



1. Physical splitting 



As explained above, metric perturbations can be decomposed according to their spin with respect to a local rotation 
of the coordinate system. This leads to density modes, vorticity modes, and gravitational waves. Gravitational (spin 
2) waves are "true" degrees of freedom of the gravitational field in the sense that they can exist even in vacuum. The 
number of gravitational wave modes is given by the dimension of the vector space spanned by symmetric, transverse, 
traceless rank 2 tensors in an n-dimensional space: this is ^(n — 2)(n -I- 1). In addition, when there is a non trivial 
matter content, there may be vorticity (or spin 1) modes arising from rotational velocity fields, which have n — 1 
independent components. Finally, there remain ^n{n + 1) — ^(n — 2)(n -I- 1) — {n — 1) — 2 possible density (spin 0) 



modes, which are usually represented by the two Bardccn potentials $ and 'J l6G, 61 
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More schematically, let us consider the metric perturbation around a locally inertial frame, written in synchronous 
gauge and in Fourier space considering the wave vector fc* = k6l : 



(1 








Vo 







2PE-2C 



ikV2 
ikVn 











ikV2 



ikV. 



-2C + Y,ht h^s 



i=3 



'23 



'■2n 



-2C - ht 







ikVn 



(5.1) 



-2C -h+J 



The quantities E and C describe the density modes (with the standard definition of the Bardeen potentials, one has 
$ = —C and \1/ = —dfE), the Vi {i = 2, . . . ,n) represent the vorticity modes, and the hj' {i = 3,. . . ,n) and /ij^^ 
{2 < j < k < n) represent the gravitational waves (when n = 3, these notations agree with the standard definition of 
the /i+ and modes). 



2. Geometrical splitting 



The three above types of perturbation generally do not evolve independently: even at linear order, they are coupled 
if the unperturbed spacetime does not possess any symmetries. However, for most cosmological models (including the 
ones considered in this paper), spacetime possesses some symmetries, being invariant under a certain group of global 
transformations. We consider the symmetry group SO(iV) with N < n, which is of course relevant when there exists 
a coordinate system in which N coordinates span a maximally symmetric space. 

When this is the case, perturbations may be decomposed into components which transform irreducibly under 
SO(7V)-rotations of the coordinate system. This leads to what we call scalar, vector and tensor perturbations which are 
perturbations whose spin with respect to SO(A^) is 0, 1 and 2 respectively. The main advantage of this decomposition is 
that the three new types of perturbation are now decoupled from each other, and hence are convenient when studying 
the evolution of cosmological perturbations. For example, consider an n-dimensional space with N coordinates 
(labelled by z, j, etc) spanning an iV-dimensional, maximally symmetric sub-space, with metric jij, and associated 
covariant derivative V^. The n — N remaining coordinates will be labelled by A, B, etc. In this case, the metric 
perturbations can be decomposed as 



5gij = -2Cjij-2VijE-2V^iEj^-2Eij, 

5gAB = E(AB), 



(5.2) 
(5.3) 
(5.4) 



where barred quantities are divergenceless A'^-vectors, and double barred quantities are divergenceless, traceless A'^- 
tensors of rank 2 (with respect to the covariant derivative Vi and metric 7ij respectively). With our definitions, it 
is clear that C, E, Ei^a)^ E(^ab) are scalars, Ei, E(^A)i are vectors, and Eij are tensors under SO(A'') rotations. The 
perturbed metric components can then be written as 



^goii3 = 



/o 









2k'^E - 2C 



-2ikE2 

ikE3 



ikE, 



N 



ikE, 



(A) 











-2ikEo 



-2ikE^ 



-2C 



N 

i=3 
'23 



ht 



''23 



(A) 2 



-2C 



ht 



E, 



(A) 3 







-2ikE 



N 



'■2N 



-2C - h+ 



E, 



{A)N 



\ 



ikE, 



E, 



(A)2 



E, 



{A)3 



E, 



E{AB) J 



(5.5) 



with the hi , hf^ describing Eij . Obviously, one has 

• 2 + [n — N) + \{n — N){n — N + 1) scalar degrees of freedom. 
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• {N — l)(n — + 1) vector degrees of freedom and, 

• ^{N — 2)(A^ + 1) tensor degrees of freedom. 

By definition, the tensor components are spin 2 quantities and represent gravitational waves. It is clear that when 
N ^ n, not all the gravitational waves are tensor perturbations (with respect to SO{N)): ^{n — N){n + N — 1) of 
them are actually scalar or vector perturbations. In fact, the spin of the second decomposition can be understood as 
the projection of the spin of the first decomposition on the maximally symmetric space. Therefore, density modes are 
always scalar perturbations, vorticity modes can be either scalar of vector perturbations, and gravitational waves can 
be any of the three. By comparing Eqns (5.1) and (|5.5|), it is clear that: 



• the 2 + (n — N) + ^{n — N){n — + 1) scalars decompose as 2 density modes, n ~ N vorticity modes, and 
i(n — N){n — + 1) gravitational waves, 

• the (N — l)(n — iV + 1) vectors represent (TV — 1) vorticity modes and {N — l)(n — N) gravitational waves, 

• the i(iV - 2)(iV + 1) tensors ah represent gravitational waves. 
For our purpose (n = + 1 = 4) , this reduces to 

• 4 scalar degrees of freedom which split into the 2 density modes, 1 vorticity mode, and 1 gravitational wave, 

• 4 vector degrees of freedom which go into 2 vorticity modes and 2 gravitational waves, 

• 2 tensor degrees of freedom which all represent gravitational waves. 

As expected, we have 10 degrees of freedom 5 of which are gravitational waves. This decomposition ensures that even 
in the vacuum, the scalar and vector parts of the Einstein equation will allow non trivial solutions. These are usually 
called "graviscalars" and "graviphotons" |64). This effect represents the most striking change to the physics of 
brane cosmological perturbations as compared to that of the standard FLRW case since it can occur at arbitrary low 
energy as long as the corresponding gravitational waves exist in the bulk. 



3. The brane point of view 



The brane is, by definition, described by A^ + 1 coordinates: one timelike and the N spacelike coordinates spanning 
an A'^-dimensional maximally symmetric space. For the case of one codimension, we have A^ — n — 1. The perturbed 
induced metric of the maximally symmetric space then has ^n{n — l) independent components. An important question 
is how these perturbations can interact with the bulk perturbations. It is clear that whatever the bulk matter content, 
there are at least i(n — 2)(n + 1) = ^n(n — 1) — 1 degrees of freedom which arise from the gravitational waves in 
the bulk. Therefore, one can expect that i(n — 2)(n + 1) of the brane perturbations can interact with the bulk. We 
will see that this is indeed the case: the second Israel condition essentially states that the discontinuity of some bulk 
perturbations which can exist even in the vacuum describe the matter content of the brane. But this also suggests that 
one additional scalar degree of freedom of the brane is likely not to be directly related with the bulk perturbations. 
It happens, indeed, that this extra degree of freedom physically corresponds to the perturbation of the brane position 
in the bulk, which is independent of the gravitational waves. For example, if the bulk is pure Minkowski space, one 
can consider a fixed coordinate system (as, e.g., Newtonian gauge, which is unambiguously fixed). The position of the 
brane in this coordinate system is defined independently of the metric perturbations. This extra degree of freedom 
ensures that in any situation all the brane perturbations can be related to bulk perturbations (see also the discussion 
in Ref. ||l^). One of the aims of this paper is to make the link between these two sets of perturbations. 



B. Geometrical perturbation variables 



We now make use of maximal symmetry on M. Due to rotational invariance, we can split the perturbations into 
components which transform irreducibly under rotations, i.e., into different S0(3)-spin components, which evolve 
independently to first order perturbation theory. One could then go on and split these into irreducible components 
under translations, corresponding to the expansion in terms of eigenvectors of the Laplacian on M (which is the 
Fourier transform in the case k = 0) Following the discussion of the last paragraph, the perturbed line element 
can be generally written as 

ds^ = n^(i + 2A)drf + 2anB,d'qdx' - a'^ij.^-j + /i^ )da;Ma;^ + 2nbB±dr]dy + 2baE±,dx'dy - 6^(1 - 2E±_L)dy'^. (5.6) 
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Here, the _L indices of E±i, Ej_±_ are labels. The quantities Bi and E±i are vectors on A4 which can be respectively 
decomposed into scalar (spin 0) components B, E±, and divergenceless vector (spin 1) components B^, E(±-)i, such 
that 'j^^WiBj = 7*-'Vii?(i)j = 0. Equivalently, the tensor on Ai, hij, can be decomposed into two scalars, C and 
E, a divergenceless vector, Ei, and divergenceless, traceless, tensor (spin 2) component, Eij, such that j''^ViEj = 
^^^ViEij ^ El — 0. This decomposition is 



B, = V,B + B„ (5.7) 

E±^ = V^E^ + E^^,i, (5.8) 

h,j = 2C7y+2S„-, (5.9) 

E,j = V^,Ej-,+E,,, (5.10) 

E^ = V,E + E,. (5.11) 



The indices of these A^-quantities are raised and lowered with the metric 7jj. The symmetries of the metric ensure 
that the scalar {A, B, C, E, B±, E±, E±±), vector {Bi, Ei, E^^-^i) and tensor (i?y) quantities evolve independently. 



C. Gauge invariant metric perturbations 

Let us consider an infinitesimal coordinate transformation 

a:"^x"+^", (5.12) 



with 



r = [T,L\L^), 
V = V'L + V. 

Under this coordinate change the geometrical perturbations transform in the following way: 



A A + du{nT) + IbL^ 

B, ^ B,--L, + -V,T, 
n a 

C ^ C + nnT + HbL^, 

Ei — > Ei + Li, 

E - E ■ 

B^ ^ B,- -L^ + ?r', 



n 



b' 



a , b ^ 

E±i E±i — -Lj ViL , 

b a 

E^^ ^ Sx± - KnT - dn{bL^), 

-b+'^e) ^ ('-b + ^e)+t. 



5.13) 
5.14) 



5.15) 
5.16) 

5.17) 
5.18) 
5.19) 

5.20) 

5.21) 
5.22) 
5.23) 

5.24) 



(Recal l that ' = d/drj and that ' = d/dy. In this section we will use both this notation and the du,n defined in 
Eqns (2^,2.4): we aim to do so in such a way as to keep the equations as simple as possible.) We can therefore define 



the following four scalar and two vector perturbation variables, which are invariant under infinitesimal coordinate 
transformations, also called gauge transformations in this context: 



n 



a2 



^- = A-du{aB + —E ] + I{ aE±_ + —E' , (5.25) 



b 



* = -C + H[aB + ^E ] - H{ aEi_ + ^E' ] , (5.26) 
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a a" ■ \ . „ I a 



S = - na„ (-B + -E - bdu i-^E^ (5.27) 



h 


= E±± 


+ U [aB 






a - 




= B,+ 


-E^, 


n 






a -, 


hi 







^Ej -dni^aE^ + —E'j , (5.28) 

(5.29) 
(5.30) 

The two vector variables possess two independent components (hence four degrees of freedom). The tensor variable Eij 
is gauge invariant since there are no tensor type gauge transformations, and possesses two independent components. 
All these quantities represent a generalisation of the Newtonian gauge often used in FLRW cosmologies (we can no 

longer call it a "longitudinal gauge" , as Sgo^ ^ 0). It is completely fixed by setting (^^B + ^E^ , E\ (^^E± + ^E'^ 

to and in this case one has 

5.900 = 2n2*, (5.31) 

6goi — anT,i, (5.32) 

Sg,j = 2a2($7y -5,,), (5.33) 

6go4 = nbYj, (5.34) 

Sgi4 = bahi, (5.35) 

Sgi4 = 2b^h. (5.36) 

This gauge is perfectly well-suited for describing the bulk perturbation without a brane. In the presence of a brane, 
however, things are more complicated since some of these quantities involve first or second derivatives with respect to 
the fifth dimension, and hence they are not always regular at the brane position (see Section VI C^ ). For this reason, 
other gauge choices are often preferred, but not essential^. 



D. Perturbed stress-energy tensor 



We now perturb the unit vectors J7" and iV", defined in Eqns (2.S-2.1C), which are the timelike and spacelike eigen- 
vectors normal to the maximally symmetric 3-spaces. It follows from the normalization conditions, Eqns (2.9, 2.10 ), 
that each vector has only four independent components. Furthermore, as U" and are eigenvectors of a symmet- 
ric tensor, they are normal to each other, Nf^U^ = 0. Hence at perturbed order there are only seven independent 
components which we denote by Wq, /q, w. They are defined by 



5U" 



-jiPw -A~ f3B^), -vl l-fiw + l3E^i_) 
n a b 

-7(-w + PA + B^), -/^, -l^ipw - 
n a 



E, 



(5.37) 
(5.38) 



Neglecting the metric perturbations, the quantity w represents the perturbation of the Lorentz boost /?, 



h 7 



(5.39) 



As usual, we will decompose Wq, /q, into scalar and vector components, 



(5.40) 
(5.41) 



^ In any case, there is no particular reason why the brane and bulk metric perturbations should be the same as the brane perturbation 
depends explicitly on the brane position, which is not a quantity that can be defined everywhere in the bulk. 
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Finally, in order to write down the stress-energy tensor, it is also useful to introduce the variables 



(5.42) 
(5.43) 



which have a decomposition into scalar and vector components similar to ( p.40y5.41| ). With these definitions, general 
perturbations of the bulk stress-energy tensor, STa/s, are 



SToo - n\Sp + 2pA), 

STo^ = -an{{p + P)v,- pB,-F{f, + E^,)), 

5Tm = -nb{SF + F{A-E^^)-pBi_), 

5T,j = a2(^F7,, +ny +2F(C7y +i;y)), 

5T,A = ha{{P-Y)U + F{v,-B,)-YEj_i), 

(5T44 = {5Y - 2YE^±_ - 2FB±_) . 

Here we have defined, according to ( 2.18| -2.21, ^.39 ): 

5p = 7'('5po+/5'^i^o) 

5Y - 7^(5^0 +/3''^Po) 

6P = 5 Pa, 

SF = f3j\Spo + 6Yo)- 



2Fw, 
2Fw. 



ip + Y)w, 



(5.44) 
(5.45) 
(5.46) 
(5.47) 
(5.48) 
(5.49) 



(5.50) 
(5.51) 
(5.52) 
(5.53) 



and we have introduced the anisotropic stress tensor liij, which again may be decomposed into a scalar, (divergence- 
less) vector, and (divergenceless, traceless) tensor components according to 



n 



(5.54) 



On investigation of the behaviour of these variables under the infinitesimal coordinate transformations ( 5.12| ) (see 
Appendix ^ , we find the following scalar gauge invariant variables 

(5.55) 
(5.56) 
(5.57) 
(5.58) 

(5.59) 





= V - 


a ■ 

n 


v! 


= Vi 


a - 
+ -E^, 
n 


f 


= f 




n 


= h 






= w 





4^ 



1_ 



b^' 



(5X» = SXo 



Xo{^B + ^E 
, n 



(5.60) 



where Xo is any scalar quantity (density oq, pressure Pq i etc) . The anisotropic stress tensor Ily- is gauge invariant by 
itself due to the Stewart- Walker lemma Notice that p, F, F and w are not scalars (since they depend explicitly 
on the choice of the coordinate system via the vector fields u" and n"), but we can, however, define the following 
gauge invariant variables. 



Jp« = ^''{5pl+ (5^5Y^) + 2Fw\ 
jytJ ^ -f^{5Y^ + l3^5pl) + 2Fw\ 



SP^ 
6F^ 



5Pr 



' 



Pj\5pi+6Y^) + {p + Y)wK 



(5.61) 
(5.62) 
(5.63) 
(5.64) 



As an example, the perturbed stress-energy tensor for a scalar field is given in Appendix G5 
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E. The perturbed Einstein equations 



The explicit forms of the perturbed Christoffel symbols, the perturbed Riemann, Ricci and Einstein tensors are all 
given in Appendices F4, F5, F8, where they are expressed in terms of the gauge invariant variables introduced 
above. We now write down the full perturbed bulk Einstein equations also in terms of gauge invariant variables. They 
split into seven scalar, three vector (each with two independent components), and one tensor (with two independent 
components) equations, adding up to the required 15 components of a symmetric 5x5 tensor. These equations 
are given below, where we indicate on the right hand side of each equation from which component of the Einstein 
equations they were derived and, when necessary, the term to which they are proportional. The seven scalar equations 
are 

A(2$ + h)+ 6K^ 
-3 {2H^ + 2HU) * - mduh - 3 (2H + U) du^ 
-3iHdn + AH^)h - GhdnH 
+3{dn + 4ff)5„$ 

+3 {dn + 3H + I) (HS) = K5 (V - FE) {00} , (5.65) 

-{U + 27{)-^ -{du+U-n)h-2du<^ ^ K5{(P + p)av^-Faf) {Ot} , (5.66) 

- 2K^ 

+2 {{U + 2n){du+U) + m^) ^ + 2^du{U + 2n) 
+2 ((/ + 2H){dn + + 3i/2) h + 2hdn{I + 2H) 

+ {d^+U)d^{h + 2^) 
-{U + 2H)du{'^ -h~Z<^) 

+(a„ + /)a„(*-2$) 
+(/ + 2H)a„(*-/i-3$) 

-\{dndu + dudn + Idu + Udn)^ 

- {{U + 2n)dn + (/ + 2H)du) S 
-S ((9„ + I){U + 2H) + {du +U){I + 2H)) 

-3 {Udu - Idn) $ - 6HHY. = [SP^ + fAH) {ij}oc 7,^- , (5.67) 

$ - + = Ksa^n {ij}(x Vy , (5.68) 
- 3 {{dudn + {H^ I)du + Hdn)^ + ndn-^ - Hd^h) 

-3Sa„W - Qa + 3(7^2 _ ^ {SFi + F{^ ~ h)) {04}, (5.69) 

^{du + n + 2U) E 

-{dn+I-H)-^ -{2H + I)h + 2dn'i' = {Favi + {P - Y)af) , (5.70) 

- A(2$ -^)- 6K<^ 
+m {du + 47i) ^' + G^duH 
+3 {du + m) 
+3 {Hdn) * + 3 {2H^ + 2HI) h-3 {2Hdn + /5„) $ 

-3 {du + 3H + U) (i/S) = K5 - FJ:) {44} . (5.71) 



The three vector equations are: 



'^{A + 2K)T.^ 



2 



^{dn+4.H){{dn+I-H)t,) 

-^{dn+AH){{du+U-nyh,) = K,(^{P + p){vl-t,)-F{fl + h)) m, (5.72) 
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{du + 2W + - {dn + 2H + I)h = Kgaflz {ij} , (5-73) 

+ ^idu + m){{dn + I-H)t,) 

-\{du + AH){{du+U-n)hi) = nz(F{vl-%) + {P-Y){Jl + h)) M , (5.74) 
and the tensor equation is 

^{A-2K)E,, + {du + m+U)duE,j-{d,, + 2,H + I)drMj = {ij} ■ (5.75) 



As a small aside, it is interesting to check our analysis of Section V A. We shall take for simplicity an empty, 



Minkowski bulk, so that the terms proportional to K , Ti, U, H and / vanish. Then the above equations reduce to 



A(/i + 2$) 


= -392$, 


(5.76) 


A(* - 2$) 


= -39^$, 


(5.77) 


AE 


= -6dudn<^, 


(5.78) 


{dl - dl - A) $ 


- 0, 


(5.79) 


dnhi 


— du^i, 


(5.80) 


{dl -dl-A) 


= 0, 


(5.81) 


{dl dl - A) E,, 


= 0. 


(5.82) 



In the vacuum, in addition to the usual two tensor modes, there are one scalar and two vector degrees of freedom 
which satisfy wave equations and represent the graviscalar and graviphoton (for a total of five gravitons, as expected). 
The remaining degrees of freedom can only exist if matter is present. They describe either density or vorticity modes. 

F. Perturbed conservation equations 

We now compute the perturbed energy momentum conservation equations. Even though they do not contain new 
information, they can provide useful evolution equations for the matter content of the bulk. Here we write them 
down just for the total bulk matter. The generalisation to several components is straightforward and is given in 
Appendix Written in terms of gauge invariant variables there are three scalar conservation equations, 

{du + m + 2U){5p^ - FYs) + rndP'^ + U{5Y^ - 6p^) 
+ {dn + 3H + 21) {6F^ + F(* + h)) 
+A {{P + p)av^ - Faf^) 
-3{P + p)du^-{p + Y)duh~FduY + Fdn{'i'~h-3<i>) = {0}, (5.83) 

{du + 3H + U) {{P + p)av^ - Faf^) 
+ {dn + 3H + I) {Fav^ + {P- Y)af^) 

+SPH'^{A + 3K)a^U+{P + p)'i> + {Y-P)h + FJ: = {i} , (5.84) 

{du + 3n + 2U) {6F^ - F{^ + h)~{p + Y)Y) 
+ {dn + 3H + 2I){6Y^ - FE) - 3H6P^ - I{6Y^ - Sp^) 

+A {Fav^ + {P- Y)af^) 
+Fdu{'f-h-3<i>) + 3{P-Y)dn^ + {p + Y)d„'^ + Fd„Y= {4}, (5.85) 

and one vector conservation equation, 

{du + m + u) ((p + p){vl - %) - F{fl + h)) 
+ (5„ + m + r) (f{vI - E,) + {p- Y){Jl + hS) 

+ ]^{A + 2K)a% =0 {i} . (5.86) 
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Finally, in order to close the system, we must specify an equation of state for SP^, 6Y'^, and Ily, as functions 
either of Sp^ or of some other non dynamical variables (such as the entropy). For example, all these quantities vanish 
for a bulk con taining; non relativistic matter. For a scalar f ield, most of them are also set to zero, as is discussed in 



Appendix |G5| . We can interpret the three scalar equations ( 5.83 - 5.85| ) as the conservation equations for Sp\ v\ and 

SFK 

Notice that with these Bianchi or conservation equations, three scalar and one vector Einstein equations are redun- 
dant and can be dropped. Formally, the seven scalar Einstein equations can be split into four dynamical equations for 
the four scalar metric perturbations \&, S and h, and three constraint equations. It happens, however, that with 
our choice of variables, the splitting is not completely straightforward. For example, the {Oa} Einstein equations are 
the constraint equations for the metric components g^, 3^4, 344, the rest being the evolution equations. However, in 
terms of our gauge invariant variables, Eq. (5.68) is obviously a constraint equation. This is because ^ involves first 
and second time derivatives of the metric perturbation E. Equivalently, the three vector Einstein equations can be 
split into one constraint equation and two dynamical equations for the variables and hi. 



VI. BULK PERTURBATION WITH A BRANE 



In the previous section we have considered the most general perturbed 5-dimensional bulk spacetime for which there 
is a perturbed maximally symmetric space orthogonal to the fifth direction. We have seen that its dynamics can be 
described in terms of four geometrical scalar perturbation variables governed by four evolution equations and three 
constraints, two geometrical vector perturbation variables governed by two evolution equations and a constraint, and 
one geometrical tensor variable governed by one tensor evolution equation. In this section we add a brane to this 
system — that is we assume, as was the case for the background, that the bulk contains a perturbed homogeneous 
and isotropic brane as a singular source. This will introduce one new geometrical degree of freedom, the brane 
displacement, whose dynamics has to be considered in order to fully describe the perturbations on the brane. 



A. Brane position and its displacement 

The perturbed brane embedding is given by 

- a" + C°(fT"), (6.1) 
X' = ff' + C(cr''), (6.2) 
= Yb + eia'^). (6.3) 

Here e is the displacement of the brane from its background position X"^ — yb, and it is a function on the brane 
worldsheet. It is a true new degree of freedom which sometimes also called "radion" jl^, On the contrary, as we 
will soon see, the perturbations in the X^ and X' directions can be set to zero without loss of generality, as they do 
not lead to any physical consequences (e.g., to a physical "deformation" of the brane) . 

It was first noticed in |22| that, when studying brane perturbations of a Randall-Sundrum background (of Ref. ]TT[), 
using the transverse and traceless gauge in Gaussian normal coordinates, the brane position is no more at constant 
y in presence of matter sources. The presence of an e 7^ can be interpreted as a bending of the brane due to the 
presence of matter or gravitational waves. The bending e can in principle be set to zero by choosing a convenient 
set of bulk coordinates such that y — yb, since by the infinitesimal coordinate transformation defined in the previous 



section, Eq. (5.13), one has 

e^e-L^. (6.4) 

This is not, however, the most general possibility. As was noted in |36| p3| , such a gauge choice will also fix the 
gauge of some other perturbation variables. In fact, if we choose L-^ such that e is zero, + frE'^ is fixed (see 



Eq. (5.24)). In a gauge invariant approach one must keep e arbitrary. 



We now can define the perturbed vector orthogonal to the brane, _L" + S ± ". One easily obtains (See Appendix p42| ) 

5±„ = {-be, -6V,e, -bE^^) . (6.5) 

Note that no vector perturbations enter in the above formula. This is just a consequence of Frobenius theorem [ p6| . 
Also, the fact that the perturbations C° do not enter in the above expression illustrates that they do not correspond 
to any physical deformation of the brane (see further comments below). 
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B. Induced metric and first fundamental form 



In this subsection we calculate the perturbed first fundamental form which will be used in the perturbed first Israel 
junction condition in the following subsections. 

We shall first look at the perturbation of the induced metric Xab{o'°')- In doing so, it i s im po rtant to recall that the 
brane embedding in the unperturbed and perturbed bulks (respectively given by (3.1-3.3), (6.1-6.3)) are different. 
Therefore, each brane variable has two contribution to its perturbation: one coming from the perturbation of the 
variable and a second contribution due to the fact that we have to evaluate it at the perturbed brane position. We 
obtain 



^00 (ct") = 2n2 (a + C° + InC° + Ibe) , 

<5xo^(^°) = an (S, - -6 + -V,C°) , 
— \ n a / 



(6.6) 
(6.7) 
(6.8) 



where, to first order in perturbation theory, the right hand side of these equations are evaluated at 2;" = X" 
(corresponding to the unperturbed embedding). As usual, we can decompose the perturbation into its scalar and 
vector parts using the metric 7^ evaluated at the unperturbed brane position: Q = V^C -I- Q- This is possible since, 
for perturbations, the time derivative 9^ and spatial derivatives di are equivalent to derivation with respect to cr*^ and 
ct' respectively. 

A word of caution is in order here. Recall that the quantities A, C, I, H, etc., are defined in the bulk. However, due 
to the presence of the brane, they may be (and in fact, are — see below) discontinuous at the brane position. For this 
reason, they may only be evaluated on each side of the brane. Therefore we must check that the above expressions 
are consistent in the sense that they have the same value when evaluated on both sides of the brane — only if that 
is the case, they can be considered well defined on the brane. As we shall soon see, this consistency in fact results 
from the first Israel condition. Anticipating this result, the above equations allow us to define in the standard way 
the brane perturbations A, B_^^ etc., 

(6.9) 
(6.10) 

(6.11) 

(6.12) 
(6.13) 
(6.14) 
(6.15) 

Using the standard 4-dimensional perturbation theory, we construct the two Bardcen potentials, as well as the brane 
vector and tensor metric perturbations. 



A 


= A + 


C" + InC" + Ibe 


B 


= 


n a 






n 


c 


= C + 


nn(° + Hbe, 


E 


= E^ 


c, 


I, 


= E,- 




4- 


= 





* EE A - a, Us 



-E 



aE, 



-E' 



'C + HlaB 



<P- H be 



aE, 



-E' 



S, = 



E,. 



-E, = S,, 



(6.16) 

(6.17) 

(6.18) 
(6.19) 



(The two first equations are equivalent to Eqns (5.21) of Ref. 
fundamental form is 



) Finally, using Eq. (3.14), the perturbed first 



^00 



2^2 (A + 
auBi, 



Ibe), 



Hbe)"fij - 2a'^E,. 



(6.20) 
(6.21) 
(6.22) 
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Sqo4 = -nb [ B± - -e ] , 



Sqa = ha ( E^i - -V^e ) , 



5q4 



0. 



(6.23) 

(6.24) 
(6.25) 



Notice that the do not appear in equations ( 6.16 - 6.25 ): this is again related to the fact that they do not represent 
physical degrees of freedom. These above expressions can also be obtained by starting from the definition 



(6.26) 



paying attention that in the perturbed and unperturbed cases, the bulk metric is not evaluated at the same position 
(2/ = Yb + e and y = yb respectively) ^ . 



C. First Israel condition for the standard 4-dimensional perturbation variables 



Using Eqns (5.1f:-5.24j6^), the coordinate transformations of the following variables are obviously continuous as 
they do not involve derivatives with respect to y: 



-B 



A + Ihe - 


-> A + Ibe + du{nT), 


(6.27) 


B, - 


^ B,--L, + -\7^T, 
n a 


(6.28) 


C + Hbe - 


C + Hbe + HnT, 


(6.29) 


- 


Ei + Li, 


(6.30) 


Eij - 


E^j, 


(6.31) 


a' - A 


fa a^ ■\ 


(6.32) 


+ —e]- 


^ -B + —E] + T. 



Given the first fundamental form ([6.20 -3.2?: ), th e first Israel conditions imply that these quantities, which are linear 
combinations of the components ( |6.20D to ( p. 23 ) of the perturbed first fundamental form Sq^p, are continuous, 



[A] + [Ib]e 

m 

[C] + [Hb]e 

m 

E 









J 





or, equivalcntly 



m 

E, 



Ihe - I \aEx 
Hbe- H (aEi 



-E' 



-E' 



= 0, 
= 0. 



(6.33) 
(6.34) 
(6.35) 
(6.36) 

(6.37) 
(6.38) 



(6.39) 

(6.40) 

(6.41) 
(6.42) 



Hence tensor perturbations are continuous and only the two scalar quantities E, h and the vector quantity hi may 
jump. The first two equations are equivalent to 



m = m = 0, 



(6.43) 
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which ensures that the brane Bardeen potentials $ and ^ are well defined. We also see that the bulk vector and 
tensor perturbation Ej and Eij may be defined on the brane where they reduce to the standard vector and tensor 
metric perturbations of a 4-dimensional spacetime with maximally symmetric spacelike hypersurfaces. Equivalently, 
the corresponding first fundamental form can be rewritten as 

^00 = 2nH + + 2q^^d,) ^^B + ^E^ , (6.44) 

Ti — ■ [ Ob Q,^ •\ • ■ 

Sqo^ = q^j^' +qoo^^{-B+^E] +q..E^, (6.45) 

— a— - \n ) -^1 

= Hh. (4 - ^')*) ' (6-46) 

which indeed reduce to 2n^^, —^qij'S^ and 2a^{^'yij~Eij), in longitudinal gauge {B = E = 0) for scalar perturbations 
and in the gauge Ei = for vector perturbations. 



D. Regularity conditions for coordinate transformations and non-standard 4-dimensional perturbation 

variables 



So far we have given the relationship between the intrinsic brane metric perturbations, the brane displacement e, 
and some of the bulk metric perturbations. Things become a little bit more involved when we consider the other bulk 
metric perturbations that appear in dgai- 

As we already noticed, the first Israel condition does not constrain E±±_ (sec Eq. ( |6.25D ). Nevertheless, since the 
transformation law for Ej_± can be written as 



E, 



E, 



UnT 



dnibL^), (6.47) 

and since all the metric coefficients must remain finite, it follows that [bL^] — (see also Appendix ^) . Furthermore, 
as the coordinate transformation — *■ x°'+^" must be invertible, we also require [L'^] = 0. Thus if b is not continuous, 
then L^(yb) = 0: 



[b]^0 



i^(y = yb) = o. 



(6.48) 



This is the only additional requirement that one must impose on the coordinate system in the vicinity of the brane. 
Note that if Z2 symmetry is ass umed, L-^ (y) = must be imposed even though b is continuous |^ . 

As mentioned in Section [II C, the only place where discontinuities or singularities are allowed is the brane position. 
When e = 0, the brane is at y = yt. However, if the brane position is perturbed, e ^ 0, the above requirement implies 



[b] e = 0, 



(6.49) 



and hence b is not allowed to jump if the brane position is perturbed. If the unperturbed metric has a zeroth order 
discontinuity in the coefficient b, the brane position must remain at e = to first order perturbation theory (see 

Figure |^) . This statement is in fact valid even in the absence of metric perturbations (see Appendi x | ^ 

Let us now consider the coordinate transformations of the variables Sqo^ and Squ given in Eqns ( |6!2^ , |6^ , 



bE 



n 

u V,e 

a 



n 

bE^, V,e 

a 



nT', 
aL'. 



The first Israel condition also states that 



[bB^] - 



[bE^ 



^2^ 



e = [bBj_] = 0, 
V,e = [bE^,] = 0. 



(6.50) 
(6.51) 

(6.52) 
(6.53) 



Therefore, in order for the transformations (6.5C,6.51) to be continuous we have to imply that a valid coordinate 
change satisfies 



[T'] 



(6.54) 
(6.55) 



Allowed 



Forbidden 
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FIG. 4: Illustration of the constraints on the brane position e implied by the unperturbed "scale factor" b. When [b] = at 
zeroth order, then the brane position e can be non zero at first order (left panel). On the contrary, if [b] 7^ at zeroth order, 
then the brane position must be zero at first order, in order the perturbation theory to be valid, so that right panel represents 
a forbidden situation. This is due to the fact that in the situation of the right panel, an infinitesimal, first order coordinate 
change could change by a large (zeroth order) amount some perturbed metric coefficients in the vicinity of the brane. (Here 
we have represented a situation corresponding to an unperturbed bulk, but of course it also holds when it is perturbed. In fact 
the perturbation of the brane position necessarily induces some metric perturbations in the bulk.) 



These conditions ensure that T and L' admit second derivatives. This standard requirement f or any va lid coordinate 
transformation is therefore preserved even in the presence of a brane. In particular (see Eqns 5.22,5.18)), this means 

that if (^^B + ^E^^ or E^' are discontinuous (which is allowed), then (^^B + f^i 



E I and E^ cannot be transformed 



identically to zero by coordinate changes. This does not prevent the quantities <i>, ^ , S, fro m being well defined. 



However, it does mean that there may not be a coordinate system in which Eqns ( 5.31 - 5.36 ) are valid (which we 

never needed to suppose). 

The first Israel condition does not require the continuity of S and hi , see Eqns ( 3.89 6.97| ) below. Finally, from 



and using the fact that [bL^] — 0, the jump 



aEi 



aEx_ + -^E' 



-E' 



is gauge invariant. Therefore the gauge invariant quantity h given in Eqn ( 5.2^ ) may contain a singular part, 

„2 



h = E±± +U[aB 



E\-dn 



aEi 



E'\)- DE, 



(6.56) 



(6.57) 



(6.58) 



which again shows that h cannot always be a component of the perturbed metric tensor in the vicinity of the brane. 
Of course, this does not invalidate the results found previously, but simply suggests that other variables may be more 
suitable to describe the metric perturbations in the vicinity of the brane. 

After these remarks on the regularity requirements of gauge transformations in a bulk-brane system, we can now 
define some further gauge invariant scalar variables in terms of which we will express the second Israel condition. 
They will also be used to write the Einstein equations for an observer on the brane when we want to compare our 
results with the usual 4-dimensional cosmological perturbation theory. 

Let us first define the gauge invariant combination 



-^1 



62 



E' 



(6.59) 



Since (f -B^ 



E' ) can be discontinuous, is defined on each side of the brane. Note that we have 



(6.60) 
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Furthermore, we set 

(be*) = T. (6.61) 

When b is continuous, T/6 can be interpreted as the "gauge invariant brane position" — that is, the position of the 
brane unambiguously defined when (^(^jE±^ + is set to zero by a suitable coordinate change (which always 

exists ifb is continuous). 

In principle, derivatives normal to the brane are not defined for brane variables. But in what follows we will also 
use dnibe"^) which we simply define as 



9„(&e«) = C/e-a„ 



aE^ + —E' 
b 



(6.62) 



In other words, the operator 9„ acts on every metric perturbation defined in the bulk but not on e (i.e., we define 
i9„e = 0). Along similar lines, one can also define d'^{be^) and 9^T. The quantity 9„(&e'*) can contain a singular term 

because of the discontinuous part of {^E±^ + ^E'^. Therefore, it is useful to define 9„T by 



9„T = a„(6e«)+DS. 



(6.63) 



This n ew quant ity can take different values on each side of the brane, so that we can define [9nX] a-nd (9nT) following 
Eqns ( FTolJalTl ). 

We may not simply continue e into the bulk as a variable e is independent of y. This is a gauge dependent 
continuation and the definitions for dn{be^) and i9^(6e'') given above would be valid only in the gauge where e is 
independent of y. From the above expressions it is also clear that the variables dn{be^) and d^{be^) and all brane 
perturbation variables which contain these derivatives, like e.g. h below, are gauge invariant only with respect to 
gauge transformations parallel to the brane. Therefore, it is important to keep in mind that these quantities are 
defined only o n (each sid e of) the brane, although the notation 9„(6e'*), 9„T may be slightly misleading. They simply 
refer to Eqns ( |6^2[|6!6^ ) . 

Using equation ( 6.59| ), one can define several gauge invariant quantities which can also only be evaluated on either 
side of the brane (that is either at y = yb + e"^ or at y = yb + e~) : 



h 



Bi -- -B + —E 



E_L_L+l( { aB + —E 



n 



{du-U)ibei), 
Ube^h~dn{be^). 



(6.64) 
(6.65) 



By comparing the last equation to Eq. ( 6.5§| ), it appears that h does not contain a singular term (and hence it is a 
quantity that has a meaning on each side of the brane). 
Using Eqns (3.3£ J6.40| ) we have 



$ = $ - HbeK 



(6.66) 
(6.67) 



The derivatives and 9„$ are defined via Eq. ( 6.62 , |6.63 ) above. The above equations will become very useful 
when considering the second Israel conditions and writing down the perturbed Weyl tensor. 

The first Israel condition states that $ and ^ are defined on the brane, and that dn^, dn^, S, and h are well- 
defined on both sides of the brane, but it does not imply their continuity. In fact, it is their discontinuity which will 
enter into the perturbed second Israel condition. Using the above definitions we have the following relations for the 
discontinuous and the continuous parts of the gauge invariant scalar perturbation variables: 



[*] = -[/]T+(/)S, 

m = [H]T-{H)E, 

[S] - m-idu-{U))E-[U]T, 

[h] = [h] + [dnT], 



(6.68) 
(6.69) 
(6.70) 
(6.71) 



and 



(*) = *-(/>T + -[/]S, 



(6.72) 
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($) = ^+{H)T-^[H]E, (6.73) 

(S) = {^) + {du-mT + ^[U]E, (6.74) 
W = {h) + {dnr). (6.75) 
If ^2 symmetry is assumed, the first of these relations reduce to 

[*] = -[I]T, (6.76) 

m = [H]T, (6.77) 

and 

W = £+ JWI, (6.78) 

($) = $-i[i?]S. (6.79) 

These relations are very important. They allow us to move freely between the bulk (non underlined) perturbation 
variables, which are defined everywhere, and the brane-related (underlined) variables, which are well defined only on 
the brane, i.e. they are either defined on the brane, like $ and or on both sides of the brane, like be^, dn^, 
S, or h, etc. 

The main difference between the brane and bulk perturbation variables is which appears in the former. The brane 
displacement, however, is not unrelated to the bulk metric perturbations: a displacement of the brane induces metric 
perturbation in the bulk (as one could have guessed by making an analogy with a charged surface in electromagnetism) . 

E. Extrinsic curvature and second Israel condition 

We define the perturbed stress-energy tensor on the brane as 

= + SP)u^Up + 2{p + P)u^^5up) - SPgaf! " + a^S^/j, (6.80) 

where 6u^ is the perturbation of the energy velocity on the brane which is given by 

Sif = (-iA,-v\-e] . (6.81) 
\ n a n J 

(The Su^ component is determined by the condition (_L^ + S_J-^){u'^ + Su'^) = 0.) The variable Jl^p is the anisotropic 
stress tensor and it is gauge invariant by itself. 



As discussed in Appendix H7, we define the gauge invariant perturbations for the energy density and the pressure 
on the brane by 

Sp^ = 5p-p(-B + ^E] , (6.82) 
— — — \ n ' 

SP^ = SP^p(-B + irE] . (6.83) 
\n ' 

Similarly we define the gauge invariant perturbation variables for the velocity on the brane, 

= V + -E, (6.84) 

n 

a -i 



v, + -E,. (6.85) 



n 



To impose the second Israel junction condition, we need to compute SlC nfi which is the difference between the 
perturbed value of ICafi at the brane position y = yb + e, and the background value of ICap. The results are given in 



Appendix EI6. The extrinsic curvature has to be compared to the perturbation of the surface stress tensor, 

Ma/3 = §TLal3 - - ^T_5q_ap, (6.86) 
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whose components are given in Appendix H7 



The perturbation of the second Israel condition, 



(6.87) 



yields four discontinuity conditions for the scalar perturbation variables S, h, ^, and the first derivatives and 



^[^-{du+U- 2W)(6e«)] = K5(£ + P)civ} 

2 



+ ih- {du + u)u ^ K5 i^sp 

— S = Ksa^n. 

In terms of the gauge invariant bulk variables these conditions read 

- AE + 3 [d„<i> - Hh + nY.] 



(6.88) 
(6.89) 

(6.90) 
(6.91) 



i [E] + (a„ - H)E 

[dn-^ + ih-idu + u)j:] 

+T{[dn] {I)~[{du+U)U]) 



dm 



3 — 



(6.92) 
(6.93) 

(6.94) 
(6.95) 



Notice that when Z2 symmetry is imposed T never appears in these equations. 

For the vector perturbation variables we obtain two discontinuity conditions for hi and the first derivative 



1 



1 



Jdn + I-H] + - [{du +U- n)h] = ^5iP + Pm - ^»), (6.96) 

- [hi\ = Kgan,. (6.97) 
Finally, there is also a discontinuity condition for the normal derivative of the tensor perturbation variable Eij : 



dnEi 



(6.98) 



Note that the Israel conditions do not give any constraint on the {a4} co mponents o f Eg. (6.87). The above constraints 
can also be found directly from the singular part of Einstein's equations ( 5.65| - 5?75| ). This is relatively straightforward 
for the vector and tensor modes, but much more involved for the scalar part, as one must rewrite the equations using 
the underlined quantities defined above, and also because one has to consider the perturbation of the covariant Dirac 
function D. However, for completeness, this has been undertaken in Appendix |lj and we have checked that both 
approaches lead to the same result. 



F. Sail equation 



As we have seen, the junction conditions are conveniently written using t he un derlined variables $, etc. In order to 
use them, we must know beK The jump of this quantity, S, is given by Eq. ( |6.91 ). As its continuous part T represents 
the brane displacement, it is natural to seek an equation describing the brane motion. As for the unperturbed case, 
such an equation is found by taking the discontinuous part of the {44} component of Einstein's equations. This yields 

3 (H) 5P} - (/) Sj^ 
1 



3P ( 5„$ - Hh + 1-LT. + -AT 



-p(5„* + //i-(a„+W)S) = 



5Y' 



(6.99) 
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where Sy'^ corresponds to the pressure perturbation along the extra dimension as measured by an observer at rest 



with respect to the brane. The relationship between SY'^ and SY is given in Eq. (7.2t) below. Equation (3.99) is the 
typical equation for the displacement of a membrane (it involves the Laplacian of the displacement T). When going 
back to the bulk (non underlined) perturbations, Eq. ( |6.99| ) becomes, as expected, a wave equation for T: 

- dlipT) - Sndu (2pT + PT) - PAT + 2KpT 



2^ 



-{P + P)T (3 {H) {H-I) + ^p{P + p) 

+T (3 {H) [Y-P] + (I) [Y + p]) = [SY^ - 3 (H) SP^ + (/) Sp* 

+p{dn'^ + Ih-{du+U)Y) 

+2du{F)E-Eidu~NH) (F) 
+i{H)^[^{P + pf + {Y-P) 



+ {I)^{-^{P + p}p+{Y + p)). (6.100) 



Note that there is nothing which guarantees a priori that the motion of the brane is stable. Even in the simplest 
case {Z2 symmetry, k = 0, no bulk perturbation contributing to the right hand side of the above equation and brane 
stress-energy tensor dominated by a constant tension term), this equation becomes 

dl + 2-d„-W^ -2-]t^0, (6.101) 
'a a J 

and the mass term becomes negative for sufficiently fast expansion rate of the brane! 

VII. THE BRANE POINT OF VIEW 

In the previous sections we have derived the bulk perturbation equations and their boundary conditions on the 
brane. This allows us in principle to solve the full system of perturbation equations in the bulk for given initial 
conditions. From these one can determine also the perturbed Weyl tensor and the second fundamental form. 

In order to make contact with 4-dimensional cosmology in this section, we want to write the perturbed version of 
the 4-dimensional Einstein equations on the brane. As for the background, this can either be done directly from the 
perturbed bulk Einstein equations (5.65 - 5775| ), or using the Gauss-Codacci equation. 



A. Projected Weyl tensor on the brane 



The full expression of the perturbed Weyl tensor SCafSjS is given in Appendix F9. Here we write only the components 
of the perturbed projected Weyl tensor, SS^g, on the spacelike direction _L" + d -L ", written in terms of the underlined 
gauge invariant variables. We have 

+ itoo + ^£oodv){y;B + ^Ey (7.1) 
+£^,vA^B+^^e\+£,^E^^ (7.2) 
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+4- {^B + ^E^ + 2f ,(,(4 - ^^$), (7.3) 

Mo4 = ^^04, (7.4) 

6£i4 = SiSqj4, (7.5) 

M44 = 0, (7.6) 



where we have set 



5Z} = '^iA + 3K)^+\A{^-h) 

+ (du + U) (a„$ + H*) - [du + 2U-H) {duh + W*) 

- {dn + 1) (dn^ -Hh + nE + i A(6e«)^ 

- {dn + 21 -H) (5„* + Ih- {du + U) S) 

-*a„ {U-H)- {hdn - {I-H), (7.7) 

= l{- {n^ + du^) + {duh + m) + {u -n)h + {H - 1)^ 

-\dn{T,-{du+U-2n){be^))), (7.8) 
= ^($-*-2/i + (if + /-2a„)(6e«)). (7.9) 

Since, by solving the bulk equations, we in principle obtain the non underlined variables, it is useful to express the 
above components also in terms of these, 

M" = -(A + 3X)4> + iA(*-/i) 
3 3 

+ {du + u) (a„$ + n^) - {du + 2u-h) {duh + u^) 

-{dn+I) {dn^-Hh + VS) 

- {dn + 21 -H) (S„f + Ih- {du+u) S) 

-^du {u-n)- {hdn - ssj (7 - H) 

+Z'e\ (7.10) 
= \ {-{H^ + du'^) + {duh + m) + {U-n)h-{\dn + I-H)Y), (7.11) 

S£^ = l($_^_2/i). (7.12) 

o 

For the vector and tensor part of the projected Weyl tensor we have defined 

SSJ = -l{A + 2K)f:i + ^{dn + H){{dn + I-H)%-{du+U-n)hi), (7.13) 

S£? = \ {{du + 2{H-U))% + {2dn + {H-I))hi), (7.14) 

= ^{{du + m-2U)du-{A-2K) + {2dn + -iH-I)dn)E^. (7.15) 

B. Perturbed Einstein equations on the brane 

With these definitions, we can now write the projected perturbed Einstein equations on the brane. They split into 
four scalar equations, 

2{A + ZK)^ 



-m{n^+dum = ^«i(E£b)E^l 



-2 {H) {3dn^ - 3Hh + 3HS + A6e») 



2 



(A + 3/s:)$) 

o 

+2 (du + SH) (H* + 5„t) 



+2*a„W = ^«^fE£b)E^b 

\ b / b 



-b 

-2(F) (a„i^ + 7A-(a„+zY)E) 

+2 (i7 + /) ( -Hh + HE+ i Afoe" 



+ ^/«5E('^^b + ^Pb+3.5f! 

B 

+r^(E(£b+a))E^^^(^" 
t-* = J«i(Ea)«^Enb 

- (iJ + 7) T 

-r^ (E(i^b+a))«^Enb 

+ |/C5a2^(nB) + (^n), 



equations, 



- i (A + 2i^) Si 



|«5 E ((^B + PB){vf » - S,) - T^bU^ * + /i,)) + (M 



3 B 



(a„ + 2H)s, = ^'^i(^Eaj«Ea 

+ (77 + 7) (hi) 



^^b 

b 



(E(^b+a))aE2' 

+|«5aE(nf) + (Mf), 
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and one tensor equation 



[du + m) duE.j - (A - 2K) E, 



{H + I)(dnE. 



m 



(7.22) 



where we have set, for the bulk matter quantities evaluated at the brane position, 



dp 


= V + 


_ 2F-eK 
n 


(7.23) 




= jpM 




(7.24) 




= 5FU 


n 


(7.25) 




= SY^ H 


-r'e« -2F-e\ 
n 


(7.26) 


a/ 




beK 


(7.27) 



These corrections follow from the fact that we have to go in a coordinate system which follows the brane. The terms 
proportional to X'e" are here because we consider the bulk matter content at y = yb + e rather than at y = yb, the 
terms proportional to come from the fact that we also perform a Lorentz boost in order to follow the brane motion, 
and the term be" in the last equation comes from the fact that the brane is bent. 

As for the unperturbed case, the continuous parts of the bulk stress-energy tensor and of the projected Weyl tensor 
appear on the right hand side of these equations, as well as the components of the continuous part of the perturbed 
extrinsic curvature. These are related through the discontinuous part of the Einstein equations to the discontinuity 
of the perturbed Weyl tensor and of the bulk perturbed matter content. The corresponding equations can be found 
in Appendix fq. 



C. Perturbed conservation equation 



The brane matter conservation equations follow from the singular part of Z)^T^" = or from the Bianchi identities. 
One obtains (see Appendix |^) 



duSp^ + m{Sp^ + SP*) 
+ (£ + p) AaiZ* - 3(P + p)du^ 
(du + Sn) {{P + p)av^) +5P" 
+ ^ (A + 3K) a^n + {P + p)* 
(du + m) ({P + pM-^r 



SF' + 



Fav^ + {P~ Y)af 



F{vl-%) + {P~Y){Jl + h) 



(7.28) 
(7.29) 

(7.30) 



Again, when there are no discontinuities in the bulk matter perturbations, one obtains the usual conservation equa- 
tions. 



VIII. CONCLUSION 



In this paper we have derived gauge invariant cosmological perturbation theory in braneworld scenarios with one 
codimension. The unperturbed background system we considered (Sections pl]-pv|) consists of a 5-dimensional bulk 
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spacetime with a maximally symmetric 3-dimensional subspace of curvature k, containing arbitrary (possibly inter- 
acting) matter with energy- momentum tensor Tq,^, and a homogeneous and isotropic 3-brane again with arbitrary 
stress energy tensor Taf3- We have not assumed Z2 symmetry across the brane. As such, our work generalises that 
of previous authors who have considered perturbation theory mainly in the Z2-symmetric case, and with specific 
bulk (and brane) matter (e.g., a bulk cosmological constant |Q or scalar field (24j). We believe that the general 
setup considered here is a necessary component of any serious attempt which may be made to tackle such important 
questions as the cosmic microwave background anisotropics in braneworlds. 

The only coordinate choice we have made is to fix the unperturbed brane to be at a given position yb in the extra 
dimension. The bulk metric is explicitly time-dependent. When the bulk contains only a cosmological constant, this 
is not the most natural coordinate system: there one would work with (static) Schwarzschild-AdSs and a dynamical 
brane 5^, |5^, However, in the case of arbitrary bulk matter and especially for the study of perturbations, we 
have found it more convenient to work in a coordinate system in which the brane is at rest. 

I n Sections we derived all the relevant background equations, ending with the brane Friedmann equation (4.8- 

4.9). As discussed in Section when Z2 symmetry is not assumed, one has additional contributions to the 4- 
dimensional Einstein tensor on the brane. In order to study these terms one has to include equations for the extrinsic 
curvature. 

In the remainder of the paper we studied perturbations of this system, setting up a completely gauge invariant 
formalism. Section ^ contains a general discussion of the classification of perturbations in an n -I- 1-dimensional space 
time, as well as the interplay between bulk and brane perturbations. An important point which we note there is the 
existence of one extra scalar degree of freedom on the brane which is not a metric perturbation (although it interacts 
with some bulk metric perturbations): this is the brane displacement. In Section VI we have derived an equation of 
motion for the gauge invariant perturbation variable describing this quantity. 

In Section ^ we introduced the perturbed 5-dimensional bulk spacetime. This led to the definition of four scalar, 
two vector and one tensor gauge invariant bulk perturbation variables given in equations ( |5^25H53o| ) . Following 
the definition of gauge invariant variables for the perturbations of the bulk matter, we were able to write down 
the perturbed bulk Einstein equations in a gauge invariant manner. The perturbed brane was then introduced 
in Section In analogy with usual 4-dimensional cosmological perturbation theory, our aim was to introduce 
two scalar gauge invariant brane perturbation variables (the Bardeen potentials), one vector and one tensor metric 
perturbation. The correct definition of these variables can only be given once the perturbed brane metric and Israel 
junction conditions are used determine the brane variables in terms of the continuous part and the jump of the bulk 
perturbations. The brane variables are defined in equations (6.16-6.19). The gauge invariant brane displacement also 



enters in these definitions. Finally the perturbed Einstein equations on the brane were derived in Section VII. As in 



the unperturbed case, they contain a contribution from the projection of the perturbed bulk Weyl tensor which in 
general have to be determined by solving the bulk equations. 

Despite the fact that we have tried to present our results as clearly as possible, the formalism presented in this paper 
is technically rather complicated. This reflects the fact that we have considered a very general scenario. The corollary 
is, however, that our results should be applicable to a whole variety of different (and possibly simpler) situations of 
interest in braneworld scenarios. In a forthcoming paper, we plan to apply this formalism to a specific model and 
solve some of the perturbation equations presented here. 
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APPENDIX 

In this Appendix, we give all the necessary formulae that were used to obtain the results presented in the text. 

Here we will consider an iV + 1 + l-dimensional bulk with, A'^-dimensional maximally symmetric, spacelike hyper- 
surfaces of constant curvature k. Therefore, here a = 0, 1, . . ., M, where M = + 1, and i = 1, 2, . . ., N. We will 
consider both the unperturbed (Section ^|-^) and perturbed cases (Section ^-Q). Both the bulk and the brane matter 
content are arbitrary as well as the global geometry of the bulk. Furthermore, we do not assume Z2 symmetry. 



APPENDIX A: SOME USEFUL FORMULAE 

1. Some ensor definitions and sign conventions 

Following the definitions of ||7^, we use the sign convention ( — h-l-), that is the signature of the metric is (H — ... — ), 
and the Riemann and Ricci tensors are respectively defined by 

R'^^^s = 9^rj^,-a,r^"^ + r^;r;,-r,:r^;, (ai) 

Rap = (A2) 



The Weyl tensor is defined by 

Caf3j5 = Ra0js — [Ra-ygfis — Rasg^j + Rf3sgaj — Rfijgas) 

N{N + " gasgp^)- (A3) 

2. Brane- related metric quantities 

The induced metric is defined by 

Xah = g^. daX'' di,X'\ (A4) 
The metric can be projected back in the bulk to give the first fundamental form 

= x^^SpA:" dqX^. (A5) 
More generally, any tensor X_^_^ defined for the brane can be projected back in the bulk using 

^ai...a,. ^^pi...p„ ... 9p„X"". (A6) 

In particular, for the brane Riemann tensor, one has 

^ Q^x" Q^xf^ Q^x-f dsX^ Rpg^^. (A7) 
One can also define the normal spacelike unit vector _Lq to the brane according to 

±^daX^^O , ±^±'^ = -1, (A8) 
and the bulk metric evaluated at the brane position can be split into 

gap = gap - J^aP, (A9) 

with 

±c.p = ±aXp- (AlO) 
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Note that Eq. (^) implies 

^^l^,x^ = -daXf^ dbX-ijc^, + r;,±p). (au) 

One can also define the extrinsic curvature according to 

/C„^=q^^i?^±/3), (A12) 

which obeys the following relations 

= £"/3, (A13) 

= 0. (A14) 

With these definition, the brane Riemann and Ricci tensors, the brane scalar curvature and the brane Einstein 
tensor can be rewritten 

R,,,a = daX^^ dhX^ dcXP diX- {R^^p, - D^^Lp)D^,^,) - i?(^±,)i^(,±p)) , (A15) 

-KJCaf) + K^l^ti0 + (A17) 

R = R + + ICp^^JC'"' - IC^, (A18) 

-igC^p + lO^K^p + \q^p [le - /Cp./C^'') + . (A19) 



APPENDIX B: BACKGROUND GEOMETRIC QUANTITIES 



1. Metric 



500 = ^2 , 5°" = 4, (Bl) 



^n^-a'l^, , g'' = -^l'\ (B2) 
5m.. = -6' , = (B3) 



2. Notations 



a ~ n 



b 



n^- , T = - , U^- (B4) 
a no 



a' ~ n' ~ b' 



H=- , [/=-. (B5) 

a n b 



3. 



Christoffel symbols 
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^'j = ^^7i,- , r o„ = ^U, (B6) 

To% = i , Tjo = m) , r^o = u, (B7) 
= i , = m'^ , = u, (bs) 

2 9 

71 - n - 

Too = -gil , r - = --^H^H,, (B9) 

r?,- = (^)r^„ (BIO) 

T^Q=i+N'H+U , r/^^=i + NH + U. (Bll) 
The superscript N means that the corresponding quantity is evaluated using the metric 'jij. 

4. Ricci tensor 

2 2 

-Roo = -{drj+U-I)U-N{d^+n-I)n+^{dy+I-U)I + N—HI, (B12) 

2 2 

= {N-l)^ijk+^-fij{dri+NH+U-i)n-^^ij{dy+NH + i-U)H, (BIS) 



i?o„ = N{-H-HH + HI + UH), (B14) 

i?MM = -{dy+I-U)I-N{dy+H-U)H+-^{dr,+U-I)U + N—raA. (B15) 

5. Scalar curvature 

i? = -Af(Af + l) + + 2iV^ - -|(a^+Z^-i)(W + iVH) + ^(5j, +/-[/)(/ + iV^). (B16) 

6. Einstein tensor 

Goo = M^IzIl f + 4k - "^hA + nM - N^iH' + H^- HU), (B17) 



2 \ P J 62 

G, = -^M(^;«^+k-2!fl^),.,+(iv-i)ic.., 

2 2 

+U-i){u + {N - l)nj + ^J^j{^y + / - C/) (/ + (iV - , (B18) 



GoM = N{-H -HH + ni + UH), (B19) 
G„M = f + - #2") _N^{n + n^- Hi) + NHi. (B20) 



Riemann tensor 



-RoMOM = b'idr, +U - n^idy + i - U)i, (B21) 

2 2 

i?oio, = a^jrj{d^+H-I)n- ^-f^jHI, (B22) 



iioiMj = a'jij{H + nH-ni-UH), (B23) 

l2 2 

ii^iMi = a%j {dy+H-U)H- -^lij-HU, (B24) 

^yfei = -a^iliklji - liiljk) I ^ + ^ ~ "^2" ) ■ ^^^^^ 
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8. Weyl tensor 



CoMOM = ^j^nH^Z, (B26) 

CMiMj = ^ b^a^Zjij, (B28) 

2 

Cijki = -j^Tj^--^{likl3i-liiljk)a^Z, (B29) 



with 



z = ^ + ^{dr, + u-i){u-n)-^{dy + i-u){i- h). (B30) 



APPENDIX C: BACKGROUND MATTER CONTENT 
1. Unit vectors 

From the fields X'"' = 77 and X'-^^ = y, one can build two unit vectors and n": 

w„ = 



Ua = (n, 0,0), 
1 



u^u'' = 1, 



-,0,0 , 



no, = (0,0,6), 

1^ 



0,0, 



-1. 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 
(C6) 
(C7) 
(C8) 



One can define the following operators: 



As in the main text, we also use 



n 

dn = -n^df, = -dy. 



H 



H 



a 



I 
I 



dun 
n 

dnU 

n 



U = 
U 



dub 
b ' 
dnb 

b ■ 



(C9) 
(CIO) 

(Cll) 
(C12) 



2. Stress-energy tensor 
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For the bulk matter, it is more convenient to introduce the unit vector Ua which represents the bulk N + 2- velocity 
of the fluid: 

C/a = (n7, 0,-6/37), (C13) 

t/^C/" = 1, (C15) 
7^(1-/32) = 1. (C16) 

Here P represents the Lorentz boost which must be performed along the y axis in order to be in the rest frame of 
the bulk matter. As usual 7 = 1 / 1 — . Due to the symmetries of spacetime, the stress energy tensor of any 
component possesses N identical eigenvalues Pq. The other eigenvalues are po (associated to the timelike eigenvector 
U") and lo (associated to the spacelike eigenvector N°'). One has 

iV« = ("-1/37,0,-17), (C17) 



n b 

iV„ = (-n/37,0,67). (C18) 

TocP = {Po+Po)UaU0-{Po-Yo)N^Ni3-Pogai3, (C19) 

= {P + p)uaUi3-{P-Y)nan0-Pga0-2Fu(^^ni3y (C20) 



Too = nV(Po + /3'l^o)=nV , T^'' = ^P, (C21) 

Tij=a^Pojij=a^Pjij , T'^ = ^Pf^, (C22) 

Tom = -nbpj\Yo + po) = -nbF , T°'' = \f, (C23) 

no 

TMM = b^7\Yo+P^Po) = b''Y , T'^^ = ^Y. (C24) 



P = 7'(Po+/32ro), (C25) 

Y = 72(^0 +/^Vo), (C26) 

P = Po, (C27) 

= /37'(po + l^o), (C28) 



/3 F 



1 + /32 p + y 



(C29) 



Po 



Yo 



p-P^Y 
I + /32 

y - ,/3V 

1 + /3- 



2 



Po = -P. (C33) 
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For any species, 



3. Einstein equations 

{W + K- H^) + NUU - N{dnH + H^) = K^+,p, (C34) 

-{du+Vl)iU + {N-l)H) + {dn + I){I+{N-l)H) = K^+,P, (C35) 

iV(a„iJ + HH - HI) ^ N{drJ-i + nH ~ HU) = k^+^F, (C36) 

_ N{N-1) ^^2 ^ ^ _ ^2^ _ jv(a„W + 7^2) + NHI = K^+^Y. (C37) 

4. Conservation equations 



= r^'[/; - i?^'7v;, (C39) 

Qi = invi4+l3fDl,),0,-bjf{Df+PfTlj)^{nTj,0,-bDj), (C40) 

T/ = 7(r^'+/3i?iJ), (C41) 

= ^{Dl+prl), (C42) 

= 7(r/-/3i5/), (C43) 

DjJ = jiDf-^Tf), (C44) 

lr,,0,il?,'), (C45 







^Fj = ^D/ =0, (C46) 
/ / 

duPf + NHiPf + pf) +UiYf + pf) + {dn + NH + 2I)Ff = Tj, (C47) 

{du + Nn + 2U)Ff+dnYf + I{Yf+pf)+NH{Yf-Pf) = Df. (C48) 



APPENDIX D: BACKGROUND BRANE-RELATED QUANTITIES 

1. Brane position 

In general, one has the brane position X°' as a function of A?' + 1 variables a". We choose 





= a° 


X' 






= Jh 



2. Induced metric 

One first builds the unit vector orthogonal to the brane: 



(Dl) 
(D2) 
(D3) 



dXi^ 

±„ = (0,0,6), (D5) 
±" = (o,0,-^^ . (D6) 
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The components of this vector have the same functional form but possibly different numerical values when evaluated 
at y = and y = yb~- Then the induced metric is given by: 

= ga(3+LaLl3, (D7) 
ga/.i" = 0. (D8) 

goo = n", (D9) 



3. First Israel conditions 



For any quantity /, we define 



/=[/](%-yb)-i) + (/), (Dii) 



where [/] is the discontinuity of /, (/) is the continuous part of / and is the Heaviside function. The first Israel 
condition states that (;a/3(yb^) = 9a/3(yb~) '■ 

[a\ = 0, (D12) 

[n] = 0. (D13) 

In particular, this means that b is allowed to be discontinuous at the brane position. Also, the derivatives of a and n 
with respect to y can be discontinuous. 

4. Extrinsic curvature 



!Ca0 = qHc.Df.L0), (D14) 
^a^i" = 0- (D15) 

'n? ~ 

Koo = (D16) 

IC = gi'-lC^, = q^-K^, = "^(^^ + ^^)' 
lO'^K^, = ^{P + NH^), (D19) 
£2_/c^-/C^^ = ^{{N -1)H^ + 2Hi). (D20) 

5. Stress-energy tensor 

Formally, one can take a stress-energy tensor of the above form to describe the brane content, provided that we 
have 

(D21) 

(D22) 
(D23) 
(D24) 
(D25) 

with 



Po 


= Dp, 


Po 


= DP, 


Yo 


= 0, 


1 


= 1, 





= 0, 



D = -^5(y-yb). (D26) 
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The condition 7 = 1 and P = has not necessarily to be satisfied but is a consequence of th e coordinate choice to put 
the brane at rest with respect to the coordinate system. Because of the Dirac term ( D26 ), P and p depend on 77,0;* 
only. Since the stress energy tensor of the brane is strictly zero elsewhere, its eigenvectors are not defined outside this 
hypersurface. The vector _L" appears as the analog of the vector N" as an eigenvector associated to the eigenvalue 
Yq = 0. Equivalently one can define an TV + 2-velocity which corresponds to the eigenvector associated with the 
eigenvalue po. Therefore, 



T 



Tib 
a/3 — 



(n,O,0). 



(D27) 
(D28) 
(D29) 
(D30) 



tensor 



(D31) 



1 



N - 1 



N 



-P + Eh 



7. Second Israel condition 



(D32) 
(D33) 



[£0/3] — I^N +2S_alz 



-N 



+ {N~1) 



TV- 1 



N - 



p + P 



^N + lPi 



— I^N + 2P^ 



Kn + 2{P + P)- 



(DM) 

(D35) 
(D36) 
(D37) 
(D38) 
(D39) 



8. Projected Weyl tensor 



(D40) 
(D41) 



^00 



TV- 1 
TV+ 1^ 
TV - 1 



~ TV(TV + 1)"'^^*^'- 



(D42) 
(D43) 
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APPENDIX E: BRANE POINT OF VIEW, UNPERTURBED CASE 

Unless otherwise noted, all the quantities are evaluated at the brane position. The quantity (9„) X stands for the 
continuous part of 5„X at the brane position. 

1. Friedmann equation 

Taking the continuous part of the Einstein equation at the brane position, we get 

[H +K) = —p +k«+2(pb) 

-NH {U) + N {dn) {H) + ^N{N + 1) {Hf , (El) 

+ {{du + u + {N -i)n)u) 

- {dn + I){I+{N- l)H) - ^N{N - 1) {Hf . (E2) 

As such, these equations are not yet very useful because they involve many terms which are not explicit 'brane 
variables'. 

2. New Friedmann equation 

Consider the combination (W{mm} + H{Qm}) of the Einstein equations. It yields 



{du + {N + l)n) [ \ ^ (n' + K- {Hf) - -^4^y j = -(AT - {HY + HF) . (E3) 

In the case {HF) = 8^, {Y) = 0, they can be integrated exactly and we find 

^ („. , . _ , ,,,, 

3. Friedmann equations using the Weyl tensor 

In general, the Friedmann equation can conveniently be rewritten using the Weyl tensor. One has 

N{N-1) 



+ ^«^+. E + PB - Fb) + ^ {Z) , (E5) 

B 



N-1 



_N{N_^ (iJ)' -{N-1) {H) {I - H) 

N-1 ^ N-1 

+ NiNTT)''^-^ E + PB + NY.) + ^^^^ {Z) . 

4, Relationship between {ICa^) and 



(E6) 
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P{H) 



N + 1 



{NP + p){H)-p{I) = ^—f[P + p + NY] + ^[Z] 



(E7) 
(E8) 



Conservation equation 



These can be found either by taking the singular part of (C47), or by considering the discontinuity of ( |B19| ). A 
bulk energy exchange term Fb can have a singular component so that Fb — DFg^'' + [Fb] [0{y — yb) ~ 5) + (rB)- 



d^p^ + NH{P^ + pJ = Fb, 

^ ^ B ' 



(E9) 
(ElO) 
(Ell) 



One also has the sail equation 



- wEa+^(^)E£b = Et^B] 



(E12) 



APPENDIX F: PERTURBED GEOMETRIC VARIABLES 



Metric 



Sgoi = anBi 
5gij = -a^hij 
SgoM = nbB± 
5giM = baE^^ 
6g,,„ = 2b^E^^ 



Sg 



00 



r.Oi 



1 



an 



1 



nb 
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(Fl) 
(F2) 
(F3) 
(F4) 
(F5) 
(F6) 



= 

hij = 

= 

El, = 



B, 



2Ei. 



V,B 

2Cjij -r ^J^ij , 

V(i£'j) + Eij, 
V^E + E,, 

ViE± + Enii. 



(F7) 
(F8) 
(F9) 
(FIO) 
(Fll) 



All 3- vectors indices are raised and lowered using metric jij. Vi represents its associated covariant derivative and 
= ViV*. Barred vectors are divergenceless, double barred tensors are divergenceless and traceless with respect to 
7y and Vi. 



2. Infinitesimal coordinate transformation 
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Under an infinitesimal coordinate transformation a;" ^ a;" + the perturbed part of a tensor transforms as 

ST^WS: - STl\::f; + ed,Tl\::f; - T-;;£— + 5,,r^T-;;-"..,^. (F12) 

Setting 

r = iT,L\L^), (F13) 

U = V'L + L\ (F14) 

the metric perturbations transform into 

A A + f + IT + IL-^, (F15) 

n . fi 

B, ^ B,--L, + -V^T, (F16) 
n a 

C C + nT + HL^, (F17) 
Ei ^ Ei + Li, (F18) 

Eij ^ ^i,-, (F19) 

^ B^--L^ + ^T', (F20) 
n 

^ E^i-^L'i--V,L^, (F21) 
a 

^ E^^-UT - L^' -UL^, (F22) 

^ ^S+4i^)+T, (F23) 

There is one subtlety due to the fact that b may be discontinuous. We consider the above infinitesimal coordinate 
transformation. For the {mm} component, we have 

5mm ^ 9mm - Tdr,{b^) - 2bdy{bL^). (F25) 

For the coordinate change to be valid, the metric components must remain finite, therefore one must have 

[bL^] = 0. (F26) 

If b is continuous, then L"*- can be an arbitrary (continuous) coordinate transformation, but if b is discontinuous, then 
L^{ri,x^,yb) = 0. Geometrically, this is related to the fact that the coordinate system is allowed to exhibit some 
pathologies only at the brane position, but not in its vicinity. 

3. Gauge invariant metric perturbations 

Using the transformation laws for + ^E^, E^, — (^fE± + fr-E'^, it is possible to construct the following 
gauge invariant quantities: 

* = ^ _ (a, + f) (£b + ^ii^) + / + , (F27) 

^^.C^n(--B+'^E)-H(iE,^p). (F28) 













- ia, ( 

n \ 


\ n 





h = E^^+U (^^B + - {dy + if) (^^E^ + ^E'j , (F30) 

Si = Bi + -Ei, (F31) 

n 

hi = E,^,,+ '^E'i. (F32) 
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8goo = 2nH + {2gaodr, + goo) (^^B + - f/„,, I^^Ei_ + '^E'^ , (F33) 

Sgoi = --9ij^' + 500 Vi ( -B + ^e\ + gijE^ , (F34) 

S9iJ = 9ii + "^E^ .g^ ^^^E^ + ^E'^ + 25.(.(4 - 6^,^), (F35) 

6goM = nbE + goody (^^B + - g^^d^ (^^E^ + ^E'^ , (F36) 

SgiM = --gijh' - gMM^^t (-B + ^e] + g.jE^', (F37) 

Sg^M = 2b''h + (^^B + - {2g^^dy + g'^^) (^E^ + ^E'^ . (F38) 

4. Christoffel symbols 

(5ro°o = ^--^IB^, (F39) 

5To° = ViA--nBi, (F40) 

n 

^^'j = J (-2^^7., + nhi, + ^k,^ + ^V(ii?,) + ^i^i?±7.,, (F41) 

^Fo^M = A'- hlB^, (F42) 

= l-{dy + i-H)B,-~id,+n + U)E^, + l^ViB^, (F43) 

5r°^ = ^{dy+i)B^-^{dr, + 2Ll)E^^-2^hlA, (F44) 

2 2 
T) ■ 71 ~ . ri ~ 

(5rjo = — vM--(a^ + H)i?* + — (F45) 

CI CI 00/ 

1 • . 1 r? . 

^^ij = 2^)+2a^^'^^~'^^^'^' ^^^^^ 

Srjk = li^jhl + Vkh) - V'h,,) + 7,, [^HB^ - ^HEl) , (F47) 

^^^^ = -2a^^''+^~"^^^"2a^^^ + ^"^^^'+2^^'^^' ^^^^^ 

JF^ = l/,;.' + l^(V'^^^--V,£;i), (F49) 

5TLm = --{dy + H)El + —UB' + ^V'E^^, (F50) 

2 2 

•^roo = —^{dr,+U)B^ + -^{dy + 2I)A + 2-^IE^u (F51) 

-^roT = ~{dy + i + H)B^-\^^{d,+U -n)E^,-\'^ViB^, (F52) 

^r^r = (2ifi;±±7., + + - l'^i^E^J) + ^^^^7i„ (F53) 

JFo'^ = -E^^ + -IB^, (F54) 

^F,^ = -Vii;^^ + ^^Eii, (F55) 

5F^„ = -E^^' + -UB^, (F56) 
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5r;;„ = d,,{A-E^^+NC + V^E), (FS?) 
Sr/;^ = \7,{A-E^j_ + NC + V^E), (F58) 
ST/^^ = dy{A-E^^ + NC + V^E). (F59) 

5. Ricci tensor 

Prom now on, we shall write any perturbation variables mostly with the gauge invariant quantities defined above. 
The non-gauge invariant terms will be put under the form of some factors involving components of the corresponding 

unperturbed tensor multiplied by (^^B + ^E^, E' and (^fEj_ + ^E'Y 

2 

n 

dRoo = ^V^* 

+{Nn + U)dr,^ + {dl + {2U - X)dr,)h + N{dl + {2n - i)dn)^ 

^-^{dy +NH + I- U){dy^ + 2/* + 2Ih) - -^Idy{^ + h + N^) 

--^{dr,dy + Udy + {NH + I)dr, + NHI + NUH + / + U')T, 

+(^00 + 2i?ooS^) [j^B + ^E^ - {R'oo + 2Ro^dr,) (^^E^ + ^E'^ , (F60) 

SRoi = {{N -i)n+u)Vi^ + {dr,-n+u)Vih + {N -i)dr,Vi^ 

In - - 

-^-^{dy + {N-2)H + 2m 



2 a a 

i an ,^ ^ , ~ / ^ „ z ^^■.^\ la. 



+ --(V2 + K{N - l))Si - ^R,jt^ 



+ 2^i9y + {N+ l)H - U) ({dy + i- H)f:,) - -^{dy + {N+ l)H - I) [{dr,+U- n)hi) 



+i?ooVi (-B + ^e) - floMVi (^E^ + ^E' ] + RijP, (F61) 

SRij = (7yV2 + iVVy +2if(iV-l)7y)$-V,j(* + 2$-/i) 

-,2 



(2H + 2n{Nn + u-i)^ + ndrji^ + h + iv$)) 

2 

-^7y (2H' + 2H{NH + i-U)h + Hdy{^ + h + JV$)) 

+ ^7^J (fidy + Hd^ + 2NnH + + H^Y. 
+ ^{d^ + {N- 1)H + U)V^itj) -\{dy + {N- l)H + i)\ 

+ (^^{dl + {NH+U- i)d,) - ^{dl + {NH + 7 - U)dy)j {Eij - 7i,$) - (V^ - 2K)E,j 

6RoM = NHdy^ - NHdr,h + N{dridy + {H - i)dr, + {n -U)dy)^ 

+ ^4 V^S + ^(/' + i{NH + 7 - U))i: --{tl + U{NH + U - X))S 
2 6 n 

+(i?OM + RQMdr, + Roody) (^^B + - {R'o^ + RoMdy + RmmO^) (^^E^ + ^E'^ , (F63) 

SRiM = -{{N -l)H + i)Wih-{dy-H + i)Vi^ + {N -l)dyWi^ 
+ l-{dn + {N- 2)n + 2Z^)ViS 



1 h 

2 a a 



{V^ +K{N -Rijh^ 
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+ ll{dn + + 1)H - U) ({dy + 7 - - 1^(9^ + (TV + l)n - 1) ((9^ +U - H)h) 

+R^M^i [j^B + - i?„„V, (^-^E^ + '^E'^ + Rij&', (F64) 

^(dr, + Nn+U- I)(d„h + 2Uh + 2W*) + -^UdJ^ + h- N^) 

-{dl + (2/ - tj)dy)-^ - {NH + i)dyh + N{dl + {2H - U)dy)^ 
+ ^idndy + {NH + U)dy + Idr, + NHI + NUH + 1 + U')T. 

+{Rmm + 2RMMdy) (^^B + ^E^ - + 2Ro^dy) (^^E^ + ^E'^ . (F65) 

6. Scalar curvature 



5R = ^ (V^(* -h-{N-l)^- KN{N - 1)$) 

+^ (2{ii + Nn){dr,+u -i)"^ + 2{ii + N-k)^ 

+ {dr,+U-i){dr,h + Ndr,^) 

+{Nn + u)dr,{h - + N{N + i){n^^ + ndr,^)^ 

+ ^ (2{i + NH){dy + i - U)h + 2{i' + NH')h 

+ {dy+i ^U){dyh + Ndy<^>)) 

+{NH + -h) + N{N + l){H^h - Hdy^fj 

~h> (^^^^ ^ ^ ^^^^ ^ ^ 

+NH' +U' + NH + i + N{N + 1)HH) S 
7. Einstein tensor 



JGoo = -^{V^ -KN{N -l)){h-^) 

+iV (^^V^ - {NH+U)dr, + ^{d^y + (NH - U)dy)j $ 

{^^^v{h - + - + ((A^ - + 2f{U){^ + h)^ 

(ndy + NUH + W') S 

Tl 

+2Goo{^ + h) - -GoM^ 


+ {Goo + 2Goo5r,) (^^B + ^E^j - {G'^o + 2Gom5^) (^^E^ + ^E'^ , (F67) 
SGoi = V^{{U-n){•^ + h)+^r,{h-<P) + N{H^ + ^)) 
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In . _ 

-- + (iV - 2)H + 27)ViS 

+^-(V' + K{N - l))Si - -GyS^' 
2 a a 

+G00V, (^^i? + - GoMVi + ^E'^ + GijE\ (F68) 

= (7yV2-Vy)(4' + $-(/i-$)-7V$) 

2 

+ ^7y (2(Z^ + {N- l)n){dr, +U-i) + N{N - ($ + 

2 

{{dl + {U-i)dn){h - (TV - 1)$) - ((TV - 1)W + w)a,(* + h- JV$)) 

2 

+ ^7..- (2(7 + (iV - l)H)idy + i-U) + NiN - 1)H^) [h - 

{{d'y + {i - U)dy){^ - (TV - 1)$) - {{N - 1)H + i)dy{^ + h + N ^)) 

2 

(dydr, + {{N - l)n + U)dy + ((TV - 1)H + i)dr, + NiN - VjHH^ S 
+2a^7i.- (;^(^ + *)S,(W + (iV - 1)7^) + ^(^ - $)5,(/ + (TV - l)i7)) 

2 

-^7ijS ((AT - l)H' +U' + (AT - 1)A + /) 
-N[^Ud,-pdy'^^-2G,,^ 

+^ (a^ + (iv - i)7i + hi) V(iS,.) - ^ (9^ + (TV - 1)^ + 7) V(Jij^ 

+ (^(^^ + (TVT^ + 7^ - T)d,) - (V2 - 2K) - ^(5^ + (TV77 + 7 - U)dy)^ E^, 

(^^ + - (i^^ + ^ + '^^(^(^^ - ^^^^^ 

<5C?0M = l(^V^+N-{H + H^-Hi-HU)-N^{H'+H^-Hi-HU)]i: 
2 n / 

+TV(7i:ay($ + *) - Hdnih - $) + (9^5^, - 7?^^ - 79^)$) 

2 \n / 

+(Gom + GoMdn + Goody) (^~^ + ^"^^ ~' (^om + Gomdy + GMhidrj) {^^^ + 52"-^'^ ' (F70) 

^GiM = -Vi((7-^)(^' + /i)+aj,($ + *) + TV(fl-/i + $')) 
+^^(9„ + (TV-2)W + 2Z7)ViS 

+^^(a^ + (TV + i)w - 7/) ((^j, + i-H)f:i) -I^i9r, + {N + i)H - 1) (idr,+u- n)h) 

+\hv^ + K{N-l))h-^-G,,h= 

+GoMVi (^B + '^E^ - GmmV^ [^E^ + ^E'^ + Gi,^;^', (F71) 
5Gmm = ^(V2-7^TV(TV -!))(* + $) 

-TV f ^V^ -^{NH- i)dr, + (NH + i)dy] $ 
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+N (^^Hdrji^ + $) + Hdy{^ + *) + {{N - 1)H^ + 2Hi){^ + 

-N^ (Hdr, + NHH + h)t. 

-2Gmm(* + /i) + -GomS 
n 

HGmm + 2GoMdy) (^^B + - (G'„„ + 2G^^dy) (^^E^ + ^E'^ . (F72) 

8. Riemann tensor 



SRoiOj = - (^n Vy + a^n-fijdr, + -^H-fijdyj * - 2-^iJJ7y (* + h) 

2 2 

o 71 - - ~ ~ ~ ~ _ on ~ — 

+—liAHdr, +ni + UH)^ + an{dr, + W)V(iS,) - — /V(i/i,) 

+ (a^(a^ + (27i: - J)5,) - ^idy^ {Eii - 7ii*) 



+RokOoiE- - <5f*) + i?o.ofc(4 - 'J-*), (F73) 
-\j ^ k] 

+(-RojOfeVj + RoijoVk) (^~^ + r?"^) ~ (-RoiMfeVj + RoijM^k) ^'^-^-L + ^"^'^ 
+Rujk (V'E - ^S') , (F74) 
<5i?oio« = -n^{^y + i-H)V^'i'-n^i\/^h+^nb{^,^ + 2U-■H)Vi^ 

+^an{dr, + 2n-U) ({dy + U- H)^^ - ^ba{dn + 2H-X) ((9^ + J - n)hi^ 

-RoMOM^i + ^E'^ + i?oyM (i;^' - ^S^) + i?o.oj (i^^' - , (F75) 

5iioiMj = (if/i - -m') - - inbVij - —lij-HU + —^ijHlj S 

W{dr,dy + {H- U)dy + {H- i)dr,){Eij - 7i,$) 
+lan{dy + H- /)V(,%) + ha{d^ + W - U)V (Jij) 



~{R'oiMj + FtoiMjdy + RmimjOti) y-^E± + ^E' 
+RokMj {E^ - <5f $) + i?o«Mfe (i^j - ^1 , (F76) 



4 4 



a^(Vfc, + i^7fcj) - -^^IkiHdn + -^IkiHdy ) {Eji - 7^7$) 
^ j] + [ik j7] — [/c <-> I] 



Defining 
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+Ri^>^^ (^B + ^e) - f^E^ + ^E') , (F77) 



SRijOM = V, [an{dy + 1- H)Y.j - ha{d^ +U- n)hjj - [z ^ j], (F78) 



-a^niik (^v,E + ^(a, + / - - ^{dr, +u- n)h^ 



-\j ^ k] 



+ (-RMiOfcVj + RMijO^k) ( — -B + ^-^1 ~ {RuiMk^j + RMijM^k) ( ^"^-L + 

(^V'S' - ^/i') , (F79) 



(^i^OMOM = -n'^idy - U)dy^ - 2n^{dy + 1- C/)(/4') - n^Idyh 
-h^idr, - I)dr,h - 2b'^{dn +U - I){Uh) - h^Udr,^ 
+nh (dr,{\dy + + dy{\dr, + W)) E 



n n? 



+{RomOm + ^RoMOudn) [ + ) ^ (-^omOm + '^RoMOividy) ( ^^J- + 'p'^' ) ' (F80) 



SRuOMi = -b^idr, +U- ii^Vih - fe^WVi* + ]^nh{dy + 21- H)Vi^ 



-\an{dy + 2H-U) ({dy + 1- + ^baidy + 2H-i) ((9^ +11- ii)h^ 

+RMOMoyi (^^B + '-^E^ + RjOM,. (^E^' - ^h'^ + Rmjmi (e' - ^S^) , (F81) 



SRuiMj = ( -b'^Vij + a^H^ijdy + -^n-fijdr, ) h + 2-^mi^ij{^ + h) 



ba? ~ ~ ~ ~ ~ b'^a ~ - ~ - 

-—lij{ndy +ni + UH)T. - — WV(iSj) + ba{dy + H)V(ihj) 

+ (a\dl + {2H - U)dy) - ^^Udr] {% - lij^) 



t2 



-{R'tAiMj + 2RMiMjdy) ^-E± + —E' 

+ RMkMj 

{E^ -5t) + RMiMk{E'; -5^,). (F82) 



9. Weyl tensor 





= {dr,+U-I)d^Eij 


71 ~ ~ _ 

+ -(5^+W-W)V(,E,.), 


(F83) 


^(02) 
ij 


~ ~ — 77 
= {n-U)dr,Eij + - 


(W-ZY)V(,S,.), 


(F84) 




= {dy + i-u)dyE,^ 


+ Udy + i-H)V(ih,), 


(F85) 




= {H-i)dyE,^ + ^{ 




(F86) 




= {dr,+u-i)({H- 


-U){^ + h)-dr,{h-^)), 


(F87) 
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= ^dy + i-U)(^{H-i){^ + h)-dy{^ + ^)y (F88) 

= (^dr,dy + {i - H)dr, + {u -n)dy + i - h +u' -H')i:. (fsq) 



-""'NWTTf'^ + (2$ + * - h) 



N 

N-1 ( ^ c?n^ 



Ar(Af + 1) 

2 ^^-1 P(01) ^ z,(02)\ ^ /I z,(^l) ^ P(^2) 



+Cofeoj - + Co.ofe (4 - S^^), (F90) 

+~'yikan{V^ + K{N - 1))S,- + anVjV^i^k) + a^VjEik 
-\j ^ k] 

+(CojOfcVj + CoijoVk) (^~^ + r?'^) ~ (CoiMfcVj + CoijMVfe) ^'^-^-l + '^^'^ 

+Cujk {V'E - ^E') , (F91) 

1 1 n^fo. 



liV- 1 



+-— ^ (an{dr^ + n-u)({dy + i ^ H)t,^ - baidr, + n-i) +U- H)h] 

-CoMOMVi (^^E^ + ^E'^ + Cj^oi (e^ - ^SJ) + CojOi [e^' - , (F92) 
+ ^an(aj/ + H- /)V(iSj) + ^6a(a^ + W - U)V(ihj) 

W (^i^dy + H- i)drj + (ia,, + n- u)dy^ e,, 

+CofcMj (i^f - + CoiMfe (4 - S']^), (F93) 
<5Cij« = --^a^ (^^k - ;^V27,fc^ 7^7(2$ + * - /i) 

^ a27ifc7,-, (V^ + i^iV) (2$ + * - /i) 



Ar2(iv + i) 



N{N + 1) '■" \n'^" ' b"^" ' nb 

n2 TV J' ^2 f^'^^^ji 
— [z .?'] + [i^ ^ i^] — [k-^l] 

-^lik^o {a^dy{^ + *) + - H){^ + h) - ^-^ {\dr, + U - W)E^ 

+^^7.fc (^{dr, + H-U) ({dy + i- H)f:,) - ^(a, + w - J) [{d^+u- n)hj)) 

+ ~l^kba{V^ + K{N - l))hj + baVjV^ihk) + a^V^J^^ 
-[j ^ k] 

(62x(0)+n2x(^)+n6XW) 



+ 



+((^01^01^ + 2ComOm9^) + ri^ ^ (^omOm + 2CoMOM9y) ^^-^J- + ^-^'^ 



N 



(an(ay + # - t^) ((^y + 7 - - 6a(9y + # - 7) ((0^ + - 



_liV- 1 
~2 N 

^^ \\Nllf '^ + (2$ + * - M 
-CMiM,(* + /i) - &'^(V' - 27C)f i,- 

62a2 f 1 ^(01) ^ p(02)A ^ , ^ 1 ^ ^^(^2) 
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APPENDIX G: PERTURBED MATTER CONTENT 



1. Unit vectors 



with 



= (^-^A,0,0^ , (Gl) 

Sua = {nA,aB,,bB±) , (G2) 

<5n" = (^0,O,~^E^j_y (G3) 

dria = {-nB^,-aE^i,-bE^^), (G4) 

SU" = (^liPw -A- 0B^), p{w + (3E^^)\ , (G5) 



SUc, = {n^{Pw + A),-a{v° -jBi-P'yE^i),-b^{w-B^-PE^^)), (G6) 
57V« = (^7(-^ + (3A + B^), -p{E^^ + pw)^ , (G7) 
5N^ = {-n^{w + pA),-a{f^ + MBi + ^E^i),b^{pw - E^^- I3B^)) , (G8) 



Pi ' 



(G9) 
(GIO) 

(Gil) 

(G12) 

(G13) 
(G14) 
(G15) 

Since rf + 6n°' is not orthogonal to u°' + 5u°', one does not have 5{n^U'^) = S{Pj), but rather 5{n^U^) = 5{f5"f) + 

2. Gauge transformation 





= 1, 




= ^1, 




Si S{pj) 


w 


Pi 1 


w 


6(3 5-f 


1 


P i 




= S{Pj), 






/° 


= v,/o + /7. 





^ Vo iL - -PlL , 

n b 


(G16) 


vf 


- - %U - IpiW, 

n b 


(G17) 


fo 


- fo + ^P^L+l^L', 
n b 


(G18) 


fl 


n 


(G19) 


w 


1 7 rji 1 l' T ± ^ T± 

W + —T+—L - -L , 
P7 P7 n 


(G20) 


6Xo 


^ 6Xo+XoT + X^L^, 


(G21) 



where Xq is any (A^ + 2)-scalar quantity (density po, pressure Pq, etc). Note: 7, /?, p, Y, F are not scalars as they 
are defined through the vector fields u", n" which specifically depend on a coordinate choice. 

3. Gauge invariant quantities 
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vl = vo + ^^E+^pjE', (G22) 
= vU-lEi + lhEi', (G23) 
/o = h-^PlE-pE', (G24) 



/T' = /? - ^/StB. - pE,', (G25) 
5X1 = SXo - Xo (^B + ^e) + X', {^E^ + ^E') . (G27) 



It is useful to define 



yi = jvi+0jfl (G28) 

vl ^ jv°^+pjf°i, (G29) 

ft = f3jvl+jfl (G30) 

ft ^ ^^^oB^^joj_ 

4. Stress-energy tensor 



6T^p = S{{Po + po)UM - 6{{Po - Yo)N^N0) - 6{Poga0) + T^ap- (G32) 

TTa/.C/" = 0, (G33) 

T^af.N^' = 0, (G34) 

TTOa = 0, (G35) 

ttmo = 0, (G36) 

Wij = a^Uij = (^(Vij - n + V(inj) + Uij^ , (G37) 

6Too = n''^\Spo + P^SYo + 2A{po+fi''Yo) + 2p{Yo + po)w) 
= n'^{5p* + 2p^) 



+(Too + 2Tood^) {-B+'^e)- (T^o + ^ToMd^) ( jE^ + ^E' ] , (G38) 



n 



b ^ 62- 



SToi = -aniiPo + pohv^ - [Yo - Po)/37/° - l^Po + P^Yo)B, - (po + Yo)E^i) 
= -an ((P + p)vl - + ^TooSi - ^Tom^ 

+TooVi i^B + ^e \ - ToMVi (^E^ + ^E') + TijE\ (G39) 



5Tij = a'^idPojij +2Po{Cjij + Eij)) 



9 \ / 2 



SToM = -nbj^ {(3{dYo + 6po) + {Yq + po){l + (5^)w + f3{YQ + pq){A - E^^) - {po + (3^Yo)Bi_) 
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= -nb{SF^ + F{^-h)-pE) 

+(ToM + To^dr, + Toody) (^^B + - (T^„ + Toudy + T^udr,) [^E^ + ^E'^ , (G41) 

6Ti^ = ba{j{Po-Yo)f°+Pj{Po + po)v° -j^Yo+(3^Po)E^,-Pj^{po + Yo)Bi) 



ha (fv\ + {P- + ^To„Si - '^T^^h 



+7oMVi (^-B + -^E^ - T^MVi \^-E± + -^E'j + TijE^', (G42) 

5T„M = b^^m6Yo + l3Hpo) + 2p{Yo+po)w-2{Yo + p^Po)E^±-2p{Yo+po)Bj_) 
= b"^ [SY^ -2Yh-2Fi:) 

+(T„„ + 2ToMdy) (^^B + - {T'^^ + 2T^^dy) (^E^ + ^E^ . (G43) 

^ l-(5pi -2p-^ -2FT) 

+(too - 2Tooa^ - 2r0"sj (^^B + ^i;^ - to°' + ^i;'^ , (044) 

-T'^V,- (^B + ^E^ - {T^^dr, + T^''dy)E\ (G45) 

5T'^ = (<5p«yi +ny ) 

+f^j (^^B + ^Ej - T'i' (^E^ + ^E'^ - 2T''^\Ej} - $), (G46) 



no ^ ' 



ST'^ = 1 + (P - y)/i) + ^r^l,- 



6a 

■ a _ 



+T'^Vj \^-Es_ + -^E'j - {T°''dn + T""dy)E\ (048) 
_^ _ ^j.MM/ _ 2T''"''dy - 2T^*'"dr,) (^^E± + ^E'^ . (G49) 



5. An example: a scalar field 



po = ^D^cpD^'cj^ + V 



Po = -Df,ct>D''<l)-V 



1 



TT - ^"'^ 



Pi 



1 

/ _ n2^' 
n(j)' 1 



b (j) r „2 
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with the ± sign determined by the condition Uq > 0. 



bJ' 



If. 

n b ' 



0^ 0^ no aqy 

b^ b^ n b acp 

Sp' = (P + p)(^-*j+(F-P)(^+/.)+fE + 

SY^ = (P + p)(^^-*j+(F-P)(^+/^)+fI:- 
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with 







at;" 










= 0, 


n 


= —hi. 



= /37^(^* + /.+ ^-^j+7^E, (G72) 

TTy 0. (G73) 

Here, the components /q are arbitrary as the eigenvalues Pq and Yq are degenerate. The expression chosen here is 
purely for convenience in order to simplify the following untilded components. 

(G74) 

(G75) 
(G76) 
(G77) 

6. Interaction term 

5Q° = - (jr" - r* - £)s) 

n 

SQi = lQi-{QOd^ + Q^dy)E\ (G79) 
5Q-^ = hsD^+Dh) 





{du + {N + m + U) ((P + p){v\ - %) - Fiff + hi)) 

+ {dn + {N+ l)H + I) (f{vI - %) + (P- Y){Jl + hi)) 



(G80) 



Q» = ViQ« + Qf. (G81) 
7. Conservation equations 

(a„ + Nn + 2U){5p^ - FY.) + NH6P^ +U{6Y^ - Sp^) 
+{dn + NH + 21) {6Fi + F{^ + h)) 
+A ((P + p)av^ - Faf) 
-N{P + p)du^-{p + Y)duh-Fdu^ + Fdni^ -h-N^) = 5r«+r*, (G82) 
{du + Nn + U) ((P + p)av^ - Faf) 
+ {dn + NH + I) {Fav^ + (P - F)a/«) 

+5P^ + ^^^^{A + NK)a'^U+{P + p)^ + {Y -P)h + FT, = aQ^, (G83) 

(a„ + + 2U) {bF^ - p(* + /i) - (p + y)s) 

+(a„ + NU + 2I){SY^ - FY.) - NH5P^ - I{5Y^ - 5p^) 

+A {Fav^ + (P - Y)af) 
+Fdu{^ -h- N^) + N{P -Y)dr,^+{p + Y)dn^ + FdnY = 5D^-Dh-VY, (G84) 



+ ^{A + {N -l)K)ani = Ql-Dhi-TYi. (G85) 
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8. Einstein equations 



A{{N - 1)$ + h) + KN{N - 1)$ 
-N {{{N - l)n^ + 2nU) * + Hduh + {{N - 1)W + u) 

-N {Hd„ + {N + l)i?2) h - 2NhdnH 
+N{dn + {N + l)H)dn^ 
+N {dn + NH + I) {HY,) 

i (9„ + (TV - 2)H + 21) S 

-{U + {N- 1)H)* -{du+U- n)h -{N- 

-K{N - 1){N -2)^ 
+ {2{U + {N- l)H){du+U) + N{N - 1)H'^) 

+2^du{U + {N - l)n) 
+ (2(7 +{N- l)H){dn + I) + N{N - \)H^) h 
+2hdn{I+{N-l)H) 
+ {du+U)du{h+{N -1)^) 
-{U + {N- 1)H)(9„(* -h- JV$) 

+(a„ + /)a„(*-(7v-i)$) 

+(/ + {N- l)F)a„(* -h- 7V$) 
-^(9„9„ + dudn + Idu + Udn)'^ 
- {{U + {N~ l)W)a„ + (/ + (TV - l)H)du) S 
-E ((5„ + I){U + (TV - 1)H) + (a„ +U){I+{N- 1)H)) 

-N {Udu - Idn) * - N{N - l)nHJ: 

iV$ -($ + *) + (/i - $) 

-N {{dudn + {H- I)du + ndn)^ + Hdn"^ - Hduh) 

-^AE - NT.{dun + - nU) 



K„+2 (V - FE) , 



1 



{du + {N- 2)n + m) E 



(G86) 



{{P + p)at;« - Faf) , (G87) 



(5P« + ^^^V^n ) , (G88) 



N 



(G89) 



«iv+2 (^i^" + F{^ - h)) , (G90) 



-{dn+I-H)^-{{N-l)H + I)h+{N-l)dn^ = K^+,{Fav^ + {P-Y)ap), (091) 
- (A((A^ - 1)4> - -I-) + KN{N - 1)$) 

+iv + {N + i)n^) * + 2iv*9„w 

+N {du + {N + 1)H) du^ 
+N {{Hd„) * + ((AT - 1)H^ + 2HI) h-{{N- l)Hdn + Idn) 

-N{du + Nn + U){HT,) = K^+,{6Yi -FT,) , (092) 



__(A + (7V-l)/f)E, 

-l{dn + {N + 1)H) {{dn + 1- H)Ei) 

+ 1 {9n + {N + 1)H) {{du +U- H)hi) 
{du + {N- l)n + U) % - {dn + {N- 1)H + I) hi 

hA + {N-l)K)hi 



{{P + P){vl - %) - F{ff + hi)) , (093) 
KN+2ani, (094) 



+ 2 (du + {N + 1)H) {{dn + 1- H)%) 
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-l{du + {N + l)H) {{du+U- H)hi) = (f{vI - %) + {P- Y){fl + hi)) , (G95) 
-{A-2K)Eij + {du + Nn+U)duEij-idn + NH + I)dnEij = K^+^Uij. (G96) 

APPENDIX H: PERTURBED BRANE-RELATED QUANTITIES 

1. Brane position 

X° = a° + C°K), (HI) 

X' = c7* + C(ct"), (H2) 

X" = Yb + eK). (H3) 

Under an infinitesimal coordinate change, the brane position e transforms into 

e^e-L^. (H4) 

2. Normal vector to the brane 



{U + SA^.)^ = 0, (H5) 

= ('-(B,-h),'- C-V^e - , -'-E^A , (H6) 

\n \ n J a \a Jo J 

5±a = (-n-e,-a-Vie,-bE^j) . (H7) 



Since _L" plays the same role for the brane as A''" for a bulk component, this means that formally, the quantities w, 
fi can be defined for the brane (we will note them w and f. respectively) at y = yh, 

w = -e, (H8) 

n 

f,+E^i = -V,e, (H9) 



a 



(HIO) 



or, equivalently, 



= -e", (HU) 
n 

/« = ^e«, (H12) 

/! = -hi. (H13) 



3. Induced metric 



^goo = 2n'^{A + ie), 

Sqoi = anBi, 

6qij = -2a'^{C + Ht)-iij -2a^Eij, 

6qoM = nb (^B± - ^e^ , 

SoiM = ba ^E±i - ^ViS^ , 

SQmm = 0. 

4. First Israel condition 



Or, equivalently 



With, 



m 



H 



e = 



m 
Si 



'-B+'-E 
n 



0, 
0, 
0, 
0, 
0, 

0. 



E. 



U-l\\E.^^^^E' 



0, 
0. 



one has 



^500 



n 

2 



Tl — ■ I Qj (X 

— gyS^ + 9ooV, -B + —E ) + giJ^;^ 

CL \ TXi Til 



^Qij = 2qk(^i ^ 



E: 
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therefore 



Then, 



0, 
0. 



[6]^ 0^X^(2/ = yb) = 0, 
[be] = 0. 



(H36) 
(H37) 



(H38) 
(H39) 



and 



therefore 



6^ b'^ 

bB_L e ^ bB_L e + nT', 

n n 

6^ b^ 

bE_Li Vie ^ bE±i V^e - aL^, 

a a 



[bB_ 
[bE^i] - 



6^1 

- e = [bBj] = 0, 
n J 

Vie = [bE^i] e = 0, 



[T'] = 0, 
[U'] = 0. 

5. New brane-related gauge invariant quantities 

At the brane position (or on both sides of the brane), 



e« = e- 



n n" 

2 



— 9„e = S — —d„e^, 

n n 



h = E^^+U (-B + ^e] - Ue = h- ]-dy{b€^), 
\n ) b 

V ^ A'-5,(a,+i)(^B + ^i;;)+J'e = *' + 7'e«-/(^^x + gi?')'. 



$ E 



n 



62 



(6e«) ^ T, 









'f 






b 




(0 



(H40) 
(H41) 

(H42) 

(H43) 



(H44) 
(H45) 



(H46) 
(H47) 
(H48) 
(H49) 

(H50) 

(H51) 
(H52) 

(H53) 
(H54) 

(H55) 
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m = 

[s] = [s]-1(s,-(zy))s- 

[h] = [h] + 



U 



(H56) 

(H57) 

(H58) 



& /- 4 



(E) + l(a,-(zY))T + il 



where the terms dyT_ are defined by setting e constant. 



(H59) 

(H60) 

(H61) 

(H62) 



6. Extrinsic curvature 



51C 



.00 



5]Cq 



SlCi 



SJCom 

^Klm finN 



-(*' + J/i) + 2/Coo*+£ooe 



+ (£oo + 2&oa^) {^B + ^E 



n{dr, + W)S - y + //i) + 2£oo* + (£oo + 2/Cooa^) 



£oo I ^■t'-L + 



o ^ a^ 
n 
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E' 



+la{dr, +U- H)hi -\'^{dy + H- i)% + -/CooSi 
+ICooV^ (-B + ^e] +lCi,P 

^nS - (ia^ + - U) (6e«) + ]^a{d^ +U - H)hi -^'^{dy + H - 1)% 



a - la a ■ 

-£ooSi+£ooVi i-B + ^E 
n \n 



n \n n \o Ir 



}C,jE\ 



-Vi, (be - b [^^E^ + j + aV(./.,) + + /C^e 



-7i,HS - -7,,($' - if ft) - Vi,(6e«) + aV(ift,) + y^,, 



Ko^goM; 



= 0. 



(H63) 



(H64) 



(H65) 

(H66) 

(H67) 
(H68) 
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iQ/3 



tensor 



The perturbed stress-energy tensor on the brane is given in Eqns ( |6.80| - 6^ ). Formally, for any quantity X defined 
on the brane, the quantity is a bulk scalar. The quantity D is also a scalar quantity of the bulk since one is 
allowed to consider the case -X_ = constant. Its perturbation is 



SB = DE^i_ - e{dy + C7)D. 



(H69) 



Note that this derivation is a bit formal: since E±± and h can be discontinuous, this expression and the next one are 
ill-defined, even in the case where b is continuous. But again, all the pathological terms cancel each othe r whe n one 
writes the Einstein equations, so that we consider that this is not a serious problem. Using the formula (G21), it is 
possible to build the gauge invariant counterparts of both 6D and SX : 



6X -X 



-B 



(H70) 
(H71) 



This last quantity is invariant under any infinitesimal reparametrization of the (t°. Equivalcntly, one has, us- 
ing (G22,G23), as well as the fact that 7 = 1 for the brane. 



With these definitions. 



a ■ 

-E, 
n 

a i, 
n 



(H72) 
(H73) 



Moo = 



N ~ 1 



N — 
-25oo3L+(f5oo 



Sp* + SP^ 



25ooa^) -B 



—an{P_ + p)v, 
a 



SSi,, = 



6S_MM 



H — ^00 
n 



5ooV, -B 



+5„ ( -B + ^E 
0. 



tE 



'^S_k{i{E^^^ 



a2 . 



(H74) 



(H75) 



(H76) 

(H77) 

(H78) 
(H79) 



8. Second Israel condition 



- AS + iV [a„$ - /i] + iVH [E] = (H80) 

\%'{du+U-2H){b(^)\ ^ n^UP + pWK (H81) 

[a„^ + ih] - [{du + um = [5P} + ^^^") ' (H82) 

-S = K„+.a^n, (H83) 

-\[d^ + I-H]% + ]^ [{du +U- 'H)K] - K^UEL + - S,), (H84) 

- \hi\ = Kjv+aafli, (H85) 
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(H86) 



9. Projected Weyl tensor 



As this quantity is defined on the branc, it is more convenient to express it in term of the brane-related (underhned) 
metric perturbations instead of the bulk-related (non underlined) metric perturbations as it was the case for the Weyl 
tensor. 



6Z^ = 



6£^ = 



+ {du + U) (6>„ $ + W*) - {du + 2U-H) {duh + U^) 

- {dn + 1) (^dn^ -Hh + H^+ A(6e«) 

- (5„ + 21 -H) (a„^^ + Ih- {du + U) E) 

(U-n)- [hdn - Wu) {I-H), 

iV- 1 



N 



(- (H* + 5„$) + (duh + U:^) + {U -n)h + {H 
-\dn{^-{du+U-2H) (6e«))), 



= j^{{N -2)^-^-{N -l)h + {{N -2)H + I-{N -l)dn){he'^)) 



(H87) 

(H88) 
(H89) 



Mr = -\^{^ + {N -l)K)% + ]^^^j^{dn + H){{dn+I-H)ti-{du+U-H)hi). (H90) 
= l((5„ + (iV-l)(7^-W))Si + ((iV-l)5„ + (iJ-/))/ii), (H91) 



= j^{{du + Nn-{N- 1)U) 5„ - (A - 2K) + {{N - + NH-I) 5„) , 



(H92) 



Moo = 



Oi 



5£, 



Om 



N - 1 



TV- 



-nWi5£'' - an5£'^ 



+fooV. [^Ib+^^E]+£,,E^, 

N -1 
N{N+1) 



0. 



(H93) 



(H94) 



(H95) 

(H96) 

(H97) 
(H98) 



APPENDIX I: BRANE POINT OF VIEW, PERTURBED CASE 
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1. New Einstein equations 

We rewrite the perturbed Einstein near the brane, that is near j/ = + e. We first define some new bulk matter 
content perturbations by 

J^^ = Spi + p'e^ -2F-e\ (II) 
n 

jp" = Jpf+P'e", (12) 

IF^ = (5P«+P'e«-(p + F)-e», (13) 

n 

= 5Y^ + Y'e^ -2F-e\ (14) 
n 

af = ap - be\ (15) 

5T^ = jr'+r'e" -£>-ef + F^5-6e«, (16) 
n a-' 

aQ^ = aQ^ + DbeK (17) 

(5£>" = 5D^-Ti: + D'e^ + (Y-P)—beK (IS) 

The terms proportional to e" come from the fact that we are considering the bulk perturbations at y = yb + e instead 
of y = yb- The other terms come from the fact that the brane is not at rest with respect to the bulk coordinate 
system, and arc a mere consequence of a Lorcntz boost with velocity = along the y axis. Going from the non 
underlined (bulk) metric perturbation to the underlined (brane-related) metric perturbations, the Einstein equation 
simplify a little bit to 

{N - 1) (A + NK) $ - N{N - l)n (W* + a„$) 
-NU (W* + du^) - NH [duh + U^) 
+Ah - Nhdn (H) + NYSuH 

+N{dn + {N + l)H)(^dn^-Hh + nE + j^Abe^^ = k^+, (d^p« + 

-(iV-1) {H^+du^) 
-{duh + U^) - {U-'H)h-{H -I)Y, 

+ i {dn + NH) {-^-{du+U- 2H) (6e«)) = {p{P^ + pja^ 



(19) 



■b 



+(Pb + pB)avl - F^af]^ , (110) 



- {N -2){N -1)K^ 
+{N - 1) {du + NH) (W3L + a„t) + {N- imdu-H 

+{N~l)U (H* + a„$) + WuU 
+ {du + 2U+{N~ 1)H) {duh + W£) 
-{N-l){d„ + I + NH) (9„t -Hh + H^) 
+ {dn + 2I+{N- 1)H) {dn^ + ih- {du + U) S) 

+hdn {I + {N- 1)H) - E5„ {I+{N- 1)H) = (dSpI + ^^^DV^Hb (111) 

+5Pi + ^^V^Hb) , (112) 

{N - 2)$ - 

+h + {dn + {N- 2)H + I) {be*) = k^+, {Da'^Uh + ^^'^b) , (113) 
-N {du + H) (dn^ -Hh + HT,+ A(6e«) 
-N-H {dn^ + Ih- {du + U) S) 

~A{^-{du+U- 2H) (6e«)) -N{H-I) 5„t = (If], + Pb*) , (114) 
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-l{du+ NU) {^-{du+U- 2H) (6e«)) 

- {^t + Ih-{du+U)^ 
+{N - 1) -Hh + n-E + ;^ A(6e«) 

+ (i? - /) ^ - (A + NK) (6e«) 

-(iV- 1) (A + 7Vi^)$ + A* 
+NWun + N{du + {N+ 1)H) {du^ + H^) 
+NH {dn^ + Ih- {du + U) S) 

-N {{N - l)H + I) f -Hh + n^+ A(6e«; 



.+2 (PBavl + (Pb - lB)a/L) , (115) 



(116) 



--(A+(7V-l)i^)Ei 
~ (a„ + (TV + l)if) ((5„ + 7 - if ) Si) 
+ ^ (9n + (^^ + l)-ff) {{du+U- H) hi) 



{du + {N- % - {dn + {N - l)H) hi 

+UY;i - Ihi 

^{A + {N-l)K)hi 

+1 (du + {N+ i)n) {{dn + I-H) %) 

-\{du + {N + l)H) {{du+U-H) hi) 



+ 2 (D(Pb+Pb)feL-S^) 

+{Pb + Pb)(4 i - Si) - F^{fl , + h,)\, 



Kiv+2 {0aU° + aUf 



). 



««+2 (Pb(4 i - Si) + (Pb - 1b) (/^ , + /li) 



(117) 
(118) 



(119) 



{du + NU) duE,, - (A - 2K) 
+UduEij - {dn + NH + I) [dnEi,) = ««+2 {p^ij + nj) . 



(120) 



(A sum on all the branc and bnlk species is implicitly assumed on the right hand side of these equations.) One can 
easily check that the singular part of the above equation reduces to the second Israel condition. 

2. Sail equation 

In the following a sum on all the brane and bulk species is implicitly assumed. 

N {H) 5P} - {I) 5^ 
1 



(121) 



NP { dn^ -Hh + n^+ — AT 



-p(5„* + //i-(a„+W)S) = 6Y 



dl{pX) - NHdu (2pT + PT) - APT + K{N - l)pT 
-N (^P + ^^^^pj X (25„W + {N + 1)H^) 

-(P + p)T [n {H) {H-I) + ^p{P + p)) 
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+r{N{H)[Y-P] + {I)[Y + p]) = [SY^ ~ N (H) SP^ + (I) Sp^ 

+NP{dn^ - Hh + UY) 
+p{dn-^ + Ih-{du+U)Y.) 
+i2du + Nn){{F)E)-Edu{F) 

+N{H)E[^iP + pf + {Y-P)) 

+ (/) S (-^(P + p)p +{Y + p)) . (122) 

3. Perturbed conservation equation 

They transform into 

{du + Nn + U) (ip - FS) + NHSP^ + U (W^ - FS) 
+ {dn + NH + 21) (W^ + + A)) 
+A (a(P + p)v^ - aFp) 
-N{P + p)du^ - (p + Y)duh - Fdu^ + Fdn{^-h- N^) = jr" + r*, (123) 
{du +Nn+U) (o(P + p)v^ - aPf^ 

+ (dn + NH + I) (aPf « + a(P - F)/) 

u jV — 1 It 

+^P" + — ^ (A + TVii') a^n + (P + p)* + (y - P)/i + PS = aQ\ (124) 

{du + NH + m) {iF^ - p(* + h)-{p + y)s) 

+ {dn + NH + 27) (W^ - Ps) - NHlP^ - I{W^ - Jp) 

+A {aFv^ + a{P - Y)/*) 
+Fdu{^ -h-N^) + N{P - Y)dn^ +{p + Y)dn^ + Fan's = W^- Dh, (125) 

{du + {N + 1)H + U) ((P + p){vj - SO - F{f! + h)) 
+ {dn + {N + l)H + I) (f{vI - %) + (P - Y){fl + h)) 

+ ^{\ + {N -l)K)aIli = Q\-Dhi-Tti. (126) 

For the brane components, they are obtained by considering the discontinuity of the {Cm}, {im} components of 
Einstein equation or by taking the singular part of the above equations. 

du5j^ + NH{5j^ + 5Pi) 

+(Pb + pJAa^;« - N{P^ + pjdu^ = M + Tb*, (127) 
E + Tb* = - E [^B + i^B*] , (128) 



{du + NU) [{P^ + p^)avi)+5Pi 



(A + NK) a^Hb + (Pb + pj* = aQl (129) 

E «2f> = - E [^B«^B + (^B - YB)aji] , (130) 



{du + {N + l)n) ((Pb + Pb)(^'"-^')) 



4. Einstein equations using the Weyl tensor 



{N -1) (A + NK)^ 

-N{N -i)n{n:^ + du^) 



-{N-1) {n^ + du^) 



(A*. - (AT - 2) (A + NK)§_) 
-1) {du + NH) (n^ + du^) 

+{N - 



{N - 2)$ - 



IN -I 
4 N 



-{N-1){H) (Ndn^ - NHh + NH^ + Afte") 



1 iV - 1 
4 W 

iV- 1 

2~ 
AT- 1 



N + 

{H) {Y,-{du+U-2H){be^)) 



(Ea)E^^ 



1 Af- 1 

4 iV 

-(iv - 1) (i/) (a„* + //i - (a„ + 

+(Af - 1) ((AT - 2)iJ + /) (dn^ -Hh + H^+ A6e« 



+ J^zr^>^..2 J2 {^Pb + Spi + N6Y 



N{N + 1) 



B 



17V- 1 
"""i N 

1 Af- 1 



(E(^b+a))E^l + 
4+2 (^E^b) «^E2b 



A^- 1 
iV(iV+ 1) 



4 iV 
-((A^-2)77 + /)T 



i4+2 (E(£b+pjj a'En 

E(nB) + (M" 



4 

AT- 1 



{A+{N -1)K)% = 



liV- 1 
4 AT 

AT- 1 



+ 



2 

A/'-l 



4+2 (Ea)E((^b+pj(^r-so) 



+^^««+2 E ((^B + PB)(t^f « - Si) - i^B(^^ « + hi)) + m 



+ {{N-2)H + I} (hi) 

B 

(a„ + N-H) dj,, - (A - 2if) E,, = ^ (^E ^b) E Si, 

+ ((iV-2)if + 7) (a„£i,-) 

+^--E(ft?) + (^S)- 

B 

5. Relationship between (SIC a/3) and [^q/j] 



( a„ t - A + + - AT 



+ 



N+1 



+ 



6Z^ 



-(^£b + \ ^"^ - ^fZi + HE + -AT 



+ 



N+1 
1 1 



5Z^ 



N+1 + 2 

-N{H) (P^ + pJavi 



+ 



(Pb + PB)av^ - Fbo/b 
N 1 



N -1 + 2 

{NP^+pJT = -N{{N-2)H + I)a'll^ 



+{N - l)a^ [Ub] + N- 



1 



5£' 



-p a{du+U-H)hi-\{dn + I-H) Si) 



+ 



-N{H) (£b+a)(^b.-SO 

(PB+PB)(4.-Si)-FB(/L + /i. 

N 1 



N-lK^^, ' 



-{NP^ + p^){h) = -N{{N-2)H + I)am^ 
+(Af - l)a [nf ] + iV— 



<5£ 
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-{NP^ + p^{^r,E,^ = -N{{N-2)H + I)m 



+(iV-l) 



5£\ 



(146) 



